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PRSFAGE 


Isxng cliains witli finite-range interactions do not show 

phase transitions. However, Ising chains with infinite-range 

interactions do show a phase transition, when the interaction 

potential is of a particular form, particularly interesting 

ai^ Ising chains v;ith short-range ant ifemromagnetic interaction 

and an infinite -range ferromagnetic interaction. Such systems 

have been studied by Na^le and Theumann and Hjeiye . Nagle 

studied spin-half Ising chains with nearest neighbour anti- 

ferromagnetic interaction and equivalent neighbour ferromagnetic 

interaction and found first order phase transition with 

ten^erature, both in a field and in tne absence of a field, 

for certain range of relative interaction strengths 

( J , , /j ) , Theumann and Hje^ye extended the 

' short-range' long-range-' 

range of antif erromagnetic interactions to the second neighbour 
and introduced I^ac potential for long-range interactions. 

These workers have reported a zero temperature phase diagram 
for the system with strength of interactions as parameters. 

In certain regions of their phase diagram, the Ising chain 
shows as many as four first order phase transitions, while in 
other regions the chain exhibits one , two or three phase 
transitions. 
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In Chapter I, we have presented results of a Konte Carlo 
simulation of spin-half Isin-g chains with (i) ne ares t-neighh our 
anti ferromagnetic interaction and equivalent neighbour ferro- 
magnetic interaction and (ii) nearest and next nearest neighbour 
antif erronagnetic interaction and equivalent neighbour ferro- 
magnetic interaction. ¥e have estimated thermodynamic properties 
such as internal energy, specific heat, staggered susceptibili- 
ties, subla.ttice magnetizations and correlation functions. The 
main aim of this work was to find out if different spin order- 
ings were present in the different regions of the phase diagrams 
of Nagle and Theumann and H;^ye and to study changes in spin 
ordering as a function of tenperature, for a given set of 
interaction strengths, in a given region. In (i) » we find a 
first order phase transition in the presence of a field while 
in the absence of a field, the system is in a metastable state 
and thus, does not show any phase transition. In (ii) , we 
observe first order phase transition in the presence and absence 
of a field, in one of the regicois . In the rest of the regions, 
we find blips in specific heat vs. tempea^ature curves, probably 
indicating that the transitions are of a higher order. ¥e find 
different spin arrangements, in different regions of the phase 
diagram at low temperatures, Tfe do not see any significant 
changes in short-range spin- ordering, within a region as the 
teBperature is varied, except in the case where a first order 
phase transition is observed. However, from the asymptotic 
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coirela.'fcion ftmcfeions coE^juted in one of tlie iTGgions , we observe 
that the long-range ordering of spins changes froa a four- 
sublattice structure to a three-sublattice structure at the 
transition ten^erature . 

The above studies on magnetic Ising chains were paxtly 
prompted by our interest in the origin of polytypism. It is 
well known that substances like SiC, ZnS, Cdl^ and TaS^ exhibit 
polytypism wherein different polytypjic forms differ only in 
the parameter of the unit cell (the other two remaining 
essentially constant) . This is due to different stacking 
arrangements of the layers in different specimens^ The nature 
of interactions responsible for poly type formation is not 
clearly known and the existing theories are inadequate . 

In Chapter II of this thesis, we have presented results 
of the Monte Carlo simulation of polyt3;pes based on the Ising 
chains of Chapter I. Ye have first started with the fact that 
each layer in a polytype can be classified e it her as a cubic 
layer or as a hexagonal layer. Therefore generation of a 
polytype is equivalent to the generation of a spin-half Ising 
chain. Since the properties of polytypes closely resesxible those 
of the magnetic Ising chains studied in Chapter I, we assumed 
that the interaction between layers are similar to the inter- 
action between spins in the Ising chains. Ho^'jever, to simulate 
polytypes we had to consider different growth mechanisms since 
the spin -flip mechanism used to generate magnetic Ising chain 
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states is not phyoicaily meaningful in the case of polytypes. 
Three different mechanisms to rearrange the layer sequence 
were construed viz. (i) rearrangement of a single layer (e.g., 

ABBCAB -ABA.CAB ) , (ii) transposition of two 

layers (e.g., ABCACB -ACBACB , and (iii) general 

rearrangement of two layers (e.g., ABCACB 

A3C3AB ) . It was found that the first two mechanisms 

presented ergodic problems, i.e. the number of trys required 
to move the chain from one configuration to another was 
tmmanageably large. The tldLrd mechanism did not present any 
such problem and hence appeared to be a more reasonable 
mechanism. Using this mechanism and with competing interac- 
tion of different ranges between layers, we have been succes- 
sful in generating poly types up to a repeat length of 12 layers. 


Though the interaction terms need justification, infinite 
range equivalent neighbour interaction at least is hnown in 
solid state systems, consisting of ions of different volumes. 

The method of simulation presented here can also be used for 
testing interaction models that may be put forth in the future 
for explaining polytypism. 


Some inorganic solids (e.g., IMAs, LaCoO^ and HoCoO^) 
and inorganic coup lex es in the solid state (e.g., 

[ Fe(phen) 2 (lies) ) , which contain a transition metal ion 
capable of showing a crossover in the ground state with 
change in crystal field strength, exhibit gradual or abrupt 
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spin-state transitions. These spin-state transitions are 
often ac coopanied by a change in crystal structure. Bari and 
Sivardiere put forth some models to explain such spin-state 
transitions, based on coupling of the spin states to the 
lattice, with some success. Because of considerable experi- 
mental interest in spin-state transitions, in this laboratoiy, 
we have tmdertaken a re-examination of the models of Bari 
and Sivardiere and also have investigated some new models for 
Spin-state transitions in solids. 

The one and two-sublattice models of Bari and Sivardiere 
do not exhibit a phase transition when a kinetic energy term 
is added to their Hamiltonian. This is an exact result. 
Similarly a dynamic model in which the spin-states are coupled 
to an optical mode, with the coupling capable of mixing spin 
states, also does not exhibit a phase transition, though it 
explains non-zero high-spin state population obseived at 
absolute zero. It appears from these results that kinetic 
energy pl?ys an important role in determining the existence 
of a phase transition. To include the kinetic energy term 
and at the same time obtain a free energy expansion in teimis 
of an order parameter we employed the path integral method for 
computing the partition function. The results from this 
method for the two-sublattice model predicted a phase transi- 
tion, in apparent contradiction to the exact results. To 
resolve the discrepancy we have examined the order parameter 
method closely and have concluded that the phase transition 



rGsuits Trom a tacit inclusion of the coupling of displacements 
on different sites, and is Inherent in the carder parameter erpan 
Sion. To explain the obssarved first order phase transitions 
we harv© considered coupling of the spin states to the cube of 
the lattice displacement and proceeding by the path integral 
method and order paraiteter expansion, we indeed observe first 
order phase transition in the model. ¥e have also investigated 
a two phonon model to satisfactorily explain the obsearved 
behaviour in the spin-state transition systems . In this 
model, the spin-states are coupled to a breathing mode and 
an optical mode. The coupling to the optical mode mixes the 
two spin-states, while coupling to the breathing mode is linear 
in displacement and without mixing, Preliminaary results show 
that this model can exhibit a phase transition and explain 
the obsearved jumps in spin-state populations, in addition to 
explaining non-zero population of the high-spin state at 
absolute zero. These results form the contents of Chapter III. 
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CHAPTER I 


MONTE GAEIO SIMUIATION OF ISING CHAINS 
¥ITH PHASE TRANSITIONS 

I«1 Introjiuotion 

One -dimensional systems are of considerable 
interest in the study of phase transitions as they are 
often amenable to rigorous mathematical treatment, 
Giirsey^ first showed that a one— dimensional gas of 
molecules interacting via a short-rdnge potential with 
negligible interaction across molecules (i.e. only 
nearest -neighbour interactions are of in^jortance) does 
not undergo a phase transition. The system, however, 
behaves like a perfect gas at high-temperatures and like 
a linear crystal near absolute zero*, for an attractive 
potential, though the partition function is analytic 
through out , A similar conclusion has been arrived at 
for a cS^^dimenSibnal system of spins represented by a 

*This is becamse the model is classical. In a quantum- 
Btechahical systemj zero-point motions preclude the 
possibility of orde^fing at O E. 
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Heisenberg Hamiltonian, in the nearest -neighbour approxima- 
tion. It is most straightforwaird to demonstrate the 
absence of a phase transition in a one -dimensional spin- 
■1/2 JLsing model with nearest -neighbour interaction. For 
the sake of sic^licity, let us consider the interaction to 
be ferromagnetic. At absolute zero, the chain is coTi 5 )letely 
ordered with all spins pointing in the same direction. 

The existing long-range order can be destroyed completely 
by introducing a single break. The free energy change in 
introducing a single break is, 

. AF = 2J - kgT In N (1.1.1) 

The first term is due to the increase in internal energy 
as a pair of spins gets alligned anti-parallel. The 
* second term is due to the increase in entropy since the 
■ break can be introduced at any of the N sites. In the 
thermodynamic limit since N-*- «► , we notice that a break can 
always be introduced at . any temperature above absolute 
zero. Thus, there can be no long-range order in the 
system at any non-zero ten^erature thereby ruling out a 
phase transition. 

h 

Mantroll has shown that for an Ising chain to 
tandergo a phase transition, it is essential that the 
highest eigen-value of the corresponding transfer matrix 
be degenerate. Since the transfer matrix will always have 
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ali positive eiements we can apply S'robenius''^ theonem 
(to see whether the highest eigen value is degenerate) * 

The theorem states that the highest eigen value of any 
matrix of finite oi*der with all positive elements is 
always nondegenerate . The order of the transfer matrix 
is finite for interactions of finite-range (in fact the 
order is 2^ x 2^ in the presence of a field where n 
implies that the interaction up to nth-neighbour is non 
zero). Eence one -dimensional systems cannot undergo phase 
transitions if the range of interaction is finite. Dobson^ 
calculated partition functions for one -dimensional chains 
interacting up to the Yth-neighbour using interaction 
potentials of different forms some of which in the 
infinite -range lead to a phase transition while others do 
not. Dobson’s calculations, however, failed to show either 
a phase transition or any significant difference between the 
two classes of potentials, thus rendering support to 
Montroll’s ideas. 

The partition function of a one -dimensional Ising 

model with infinite -range interactions was first ccrrcuted 

7 

by Eac for the potential of the following foimi: 

V(x) = + .for 0 X < 5. , V(x) = -o exp (-yx) for x > 6 
This form of the potential is now 3mown as Kao potential. 
Kac’s result showed no phase transition for finite non-zero 
Baker computed the partition function in the limit Y O 
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»Y 

or e . ■►I , There were two different situations* 

(i) lim ’♦’1 and lim e”^ *► 0 as N-*- •} (ii) ^ t 

and lim e -►1 asN-*-*. In case (i)* Baker observed a 

second o3cd.er phase transition while results of case (ii) 

9 

were equivalent to the Bragg -¥illiams approximation. 

Eac, Uhlenbeck and Hemmer^® determined the partition 
function exactly in the thenaodynamic limit for the Eac 
potential and showed that when the range of attractive 
force goes to infinity (with strength becoming porpor- 
tionately weaker) , a phase transition occurs . 

For an Ising chain with infinite-range interaction 
the Hamiltonian is given by, 

H = - 1 S.S., J(n)>0 (1.1.2) 

±>3 ■ ^ “ 

¥e can see that there can be no phase transition if 
defined by, 

Mo = I J(n) (1-1 -3) 

n=1 

is infinite, simply because there is an infinite energy 
gap between an ordered ground state and all other states. 

On the other hand, we have already discussed the 
case in which only finite number of «J{n) are non-zero and 
arrived at the conclusion that there can be no phase 
transition and that the system is disordered at all non- 
zero temper at tares. It has been proved that the above system 
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is disordered at all finite t®3i5>0rature3 and therefore 
there is no phase transition, provided that defined by 

J n J (n) , (1.1 *k) 

n=1 

is finite. This is intuitively clear as 2K^ is the energy 
required to have a single break in the corapletely ordered 
phase. Finite implies that breaks in ordering must 
appear with non-zero density along the chain at tempera- 
tures above zero and consequently no long-range order 

1 1 

can be expected. Kac and Thompson, after studying some 

particular examples, conjectured that the above system 

with finite M and infinite has a phase transition 
o 1 

above absolute zero. This is satisfied if and only if 

j(n) has the form n”® with 1<a jf, 2 .- The existence of a 

phase transition for the above potential for 1 ® < 2 has 

11a 

been proved by Dyson • for « = 2 it is believed that 
there is no phase transition. 

So far, we have stated that in one -dimension, a 
phase transition can occur only if the interaction poten- 
tial has an infinite-range. This is true as long as we 
restrict ourselves to a pairwise interaction potential 

or a many body interaction potential to a finite order. 

12 

Fisher has exact ly solved a one -dimensional mode l in 
which there is many body interaction to an infinite order 
though the range of the interaction potential is finite. 
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He has shown thal: if 'bhe many body interaction potential 

(s refers to the number of particles interacting) 

-2 

decays slower than s , the system exhibits a phase 
transition. However » here one can argue that the role 
of infinite -range interaction is taken over by the infinite - 
order many body interaction. 

Ising systems of interest to us are one -dimensional 
systems in which both short-range and long-range inter- 
actions are operating. A system with short-range eind 
long-range interaction was first studied by Baker^^ in 
higher dimensionality . In one -dimension, Suzuki^ ^ was 
the first to study systems in which both short- and long- 
range interactions were pre^nt . Suzuki ‘s Hamiltonian 
is given by. 


-mH J 


( 1 . 1 . 5 ) 


As we can see, in the above Hamiltonian, the short-range 
interaction is between the transverse components of the 
®P^ the long-range interaction is represented by an 


infinite -range equivalent neighbour interaction term first 
studied by Tenqjerley . Another point to note is that the 
shojrt-range and long-range interaction constants are of the 


samfe sign. Suzuki solved this model exactly and showed 
the existence of a phase transition. Anderson, Yuval and 
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TTamarirt liavs studisd 3- OH© — dimensional Ising system T^itli 

competing short-range and long-range interactions. Their 

long-range interaction texmi whs f'erromagnetic and decayed 

as r**^^ They showed that two phase transitions exist in 

the system for strong nearest-neighbour antiferromagnetic 

interaction coupled with inf ini te -range interaction. 

However, approximations in their theory break down for 

the short-range interaction values considered by them^ 

Nagle has considered the sane system with equivaleht 

-2 

neighbour interaction in place of the ^ interaction of 
Anderson et al. This system showed interesting phase 
behaviour although there was only a single transition in 
zero field. Theumaim and extended the short" . 

interaction to second-neighbour and used the Kac potential 
in the limit Y ^ 0 for the inf inite-range ferromagnetic 
term. They observed that for a certain range of values 
of the interaction strengths the system showed as many as 
four first-order phase transitions at zero temperatTxre. 

In the present study, we have carried out a Monte 
Carlo* simulation of the models of Nagle and Theumann and 
H/(ye , as a function of temperature , in the different 
regions of their phase diagrams , in order to imderstand^^^ 
the changes in spin ordering accossgjanying phase transitions . 
Tfe shall discuss these two models in some detail in the 
next section. This stiady was also proi^ted to some extent 
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by our interest in long period polytypss (discussed in 
Chapter II of the thesis) which can be treated as a 
spin-1/2 Ising chain with different spin orderings. 


1.2 Ising Models of N^le and Theumann ^d.H/ye 

Nagle^"^ studied the competing interaction laodeh to 
see whether such an Ising chain exhibited sietaHiagnetisiE, 
i*e.,to see whether there exist two critical tenqjeratures 
a lower critical temperature at which the chain becomes 
ferromagnetic from a lower temperature antiferromagnetic 
phase and a higher critical temperatm’e at which the 
f erromagnetism vanishes. His reasoning v/as based on the 
fact that Ising chains with an inf ini te -range interaction 
of the form ~ show a decreasing T as cr increases, 

w 

provided the total interaction strength is kept constant . 
This iniplies tha-t longer^ r*angs interactions are more 
efrectivs at higher temperatures* Thus, an Ising chain 
vitli long-range ferromagnetic interaction and. short- 
range antif erromagnetic interaction, at low t smperat tires , 
is expected to he in an antif erromagne tic state* As the 
tes5>eratur9 is increased, the longer range ferromagnetic 
interaction takes over and the chain becomes ferromagnetic 
from an antiferromagnetic state at a lower critical 
temperature* This ferromagnetism vanishes at an upper 
critical tenperature* 



9 


Nagle *s Hamilt oiiiasi is given by, 


N 3r ^ 

^ J /'i’^i+l “ N . \ . '^i^d 

i= I 3 ,1= I 


N 

H I Vi» 


i=1 


"^Sr ^ ^Zx- ^ 


( 1 . 2.1 ) 


The free energy expression corresponding to this Hamilto- 
nian was exactly derived by Nsigle . The above Hamiltohian 
without the long-range part was considered and the free 
energy expression corresponding to this (Nagle's Hamilto- 
nian without the long-range part) is obtained as a fimction 
of teniperature and the magnetic field, . A Legendre 
transformation ^ was effected to transform the variables 
in the free energy expression from and T to M and T.' 

The contribution to free energy from the long-range part 
was calculated to be The free energy expression 

corresponding to Nagle's Hamiltonian as a function of M 
and T is thus obtained. To locate the critical points, 
free energy was expanded in a power series in Ms 

A (M,T) = a^(T) + a^CT) + a^(T) + ...(1.2.2) 

Expression for a , a. etc. are known exactly as we know 

the free energy expression. The critical points are given 

20 

by the conditions, 

(1.2.3) 


^2 ^'^c> = ^4 ^'^c)> 
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Both the conditions are satisfied for R®3»1532 

and one finds a classical transition. ¥hen -J,_/J _ <R, 
we can define a ten^erature T* at which a2(T*)=0 but 
a 2 ^(T*) < O so at a higher temperature , the M=0 state 
becomes unstable and a first-order transition results in 
zero field. At there is a confluence point and 

the critical exponents are determined for this point. 

The model also exhibits a first order transition in a 
magnetic field. 


Stell and Hemmer studied the analogous linear gas 

with a hard core, a short-range repulsive interaction and 

an infinite -range attractive interaction - and arrived at 

22 

similar results, Nagle and Bonner studied Nagle's model 
with an additional staggered field and obtained phase 
diagrams in the space of staggered and direct fields and 
interaction strengths. We shall not discuss this model 
any further since it is not relevant to the present study . 


*1 8 

Theumann and H/ye extended the range of short- 
range interactions up to second neighbour and studied the 
following Hamiltonian; 


H = - 


N. 

J I 

1=1 

- H 


N 

I 

i=1 ^ 


N 

- K I 
i=1 


^i^i+2 


N 

I 4i(i-j)ji. 
i>d=1 ■ 

(1.2.4) 


and 


= +1 , J < O, K < O 
j) = ay e~^‘'^i, a > 0. 


(1.2.5) 
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The long-range paort in the limits 0 is equivalent to the 
long-range pairt of Nagle, The free energy of this acid el 
was exactly obtained by using the magnetic analogue of 
lebowitz -Penrose solution for systems with long-range 
and short -range interactions. Lebowitz -Penrose solution 
when applied to a nagnetic system, for non-negative ♦(j) 
in the limit Y -*■ o, gives the exact free energy, A(M,T) , 
from the expression. 

A(M,T) = GB [A (M,T) - a M^] , (1.2.6) 

s 

a = 1 I <^(i) (1.2.7) 

ijs-oo 

where CE stands for convex envelope and A (M,T) is the 
free energy of the system with the long-range part being 
zero identically. . A (M,T) is obtained by calculating the 
partition function of the corresponding Hamiltonian, along 
the lines of Dobson^ and later performing a Legendre trans- 
formation to the variables (M,T). The critical points are 
determined by solving for the isotherms. ' The isotheimis 
are given by the equation 

[ W MG [H^ - 2al4] , (1.2.8) 

'■ O'' 2h 

where MG stands for M^txwell construct^’ and is the 
isotherm of the system in the absence of long-range 
interactions. 
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The zero temperature phase diagram of Theumarm amd 

H(/ye as a function of interaction strengths is shown in 

Fig. 1 . 1 . ¥e see that in some regions the system shows 

as many as four first order transitions. This was also 

2 ‘^ 

observed by Stell and Hammer '' in the case of an eina logons 
linear gas and four different typies of ordering were 
observed at zero 1 01:5)6 nature . Theumann and H/ye extended 
the phase diagram to non-zero temperatui^s without signi- 
ficant cheinges in the regions. They further determined the 
critical temperatures of the system for some given sets of 
interactions and obseirved many phase transitions in both 
zero and non-zero fields. Another interesting feature of 
this model (of relevance to polytypes) is the coexistence 
of as many , as five phases for a given set of interactions . 

This seemingly violates G-ibb’s phase rule as the system 
in question is a single con5)onent system. However, since 
the variables in the phase diagram are not strictly thermo- 
dynamic, there is really no violation of the phase rule. 

Because of our interest in polytypism and its relation 
to Ising chains (for details., see Chapter II) > we have carried 
out a Monte Carlo simulation on these Ising chains. ¥e 
have attempted to find out the '^changes in the ordering of 
spins , as a function of temperature in each region of the 
phase diagram of Theumann and H;( 5 ye. Monte Carlo simulation 
was decided upon since this method in addition to giving * 
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tliermod.yn£anic properties of the systea, allows us to taike 
a direct look at any of the systems of the canonical 
ensemble, A brief introduction to the Monte Carlo method 
and its application to estimation of thermodynamic proper- 
ties is given in the next section. This is followed by 
a description of the programs and their testing, employed 
in the present study. 

1.3 Monte Carlo Method^and_ its Applications to Thermo- 
dynamics 

Monte Carlo method ~ involves the use of random 
saa^ling to approximate the solution of maths mat ic a 1 or 
physical problems. Mathematical problems solved by Monte 
Carlo method are mostly confined to the evaluation of 
multi -dimensional integrals. Other mathematical problems 
solved by Monte Carlo method {although the method is much 
inferior compared to the standard numerical methods) are, 
determination of eigen values and eigen functions of 
matrices, solutions of differential and integral equations, 
inversion of matrices and such others, Monte Carlo method 
is most useful in the evaluation of multi-dimensional 
integrals, where the conventional numerical methods become 
unmanageable as the dimensionality of the integral becomes 
large. The error in the Monte Cfsuclo method applied to 
mathematical problems is proportional to o where o 
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is the sasipliiig variaaice of the estimator and N is the 
sa35)le size. Thus, error in Monte Carlo estimates can be 
reduced by either increasing the sample size or by 
decreasing the variance. The former method of increasing 
accuracy is often not practical as a reduction in error by 
ten-fold demands a hundred -fold increase in sample size or 
correspondingly a htindred-fold increase in coE^uter time. 
Efficient Monte Carlo methods should therefore depend more 
on reducing the variance of the estimator. In recent times, 
a combination of fast computers and good variance reduction 
techniques have made accurate estimation of multi-dimensional 
integrals by Monte Carlo method, plausible. 

Several other mathematical problems mentioned earlier 
have been solved by devising ingenious stochastic processes 
such that the distribution of the so devised stochastic 
process or some set of its parameters satisfy the mathemati- 
cal equation in question. One of the earliest problems 
solved -was the Buffon needle problem for estimating it - 
the value of ir was obtained by determining the fraction 
of times a needle of given length, when dropped randomly 
between a pair of parallel lines separated by a known 
distance, intersected one of the lines. 

Direct Monte Carlo simulation has been used to, stiady; 
(i) physical problem of a cos35>le3 nature where no analytioal 
formulation is practicable, and (ii) to gain insight into 
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the physical processes taking place in a systesa. Direct 
sisiulation involves playing a game of chance in the system, 
simulated say on a computer, with the rules of the game 
being ouch that the actual deterministic and random features 
of the physical process are step by step exactly imitated 
by the game. Relevant parameters of the physical system 
can be obtained fairly accurately when the game is played 
for a sufficient length of time. Large number of problems 
from diverse fields are being solved by direct Monte Carlo 
simulation examples being neutron transport problems and 
nuclear reactor criticality, ^ quality and reliability of 
con5)l6x systems, mass servicing systems, percolation 

26 

problems, chain conformation, molecular weight distribu- 

te 

tions and polymerization rates of polymers apd so on. 

Any Monte Carlo calculation involves the generation 

of random numbers belonging to a specified distribution. 

The random numbers so generated should satisfy appropriate 

statistical tests such as frequency tests, moments tests, 

chi-squared test etc. This is important since errors in 

the estimates can be reduced by using a well distributed 

set of random numbers . Furthermore , the time taken for 

each random number generation must be as small as possible 

and the random number sequence so generated must have long 

periodicity. There are several algorithms for fast 

31 32 ■ 

generation of well distributed random numbers * The 
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random number generation algorithm used in the present 
study as veil as the results or statistical tests carried 
out on the Tirst thirty thousand random numbers generated 
by using this algorithm are given in the Appendix at the 
end of this chapter. 

Before considering the applications of Monte Carlo 

method to estimate ther!nod 3 mamic properties, a brief review 

of Markov chains is in order. An nth order Markov 

chain is one in which the probability of outcome of a 

particular trial depends on the previous n outcomes and 

none else (for n=:0, we have independent outcomes for each 

trxal and the process is known as Bernoulli trials) . ¥e 

shall be conceimed only with first order Markov chains or 

simple Markov chains where probability of outcome of a 

particular trial is dependent on the outcome of the 

previous trial only. All possible outcomes of a Markov 

chain are said to form a state space. A Markov chain with 

a finite state space is known as a finite Markov chain. 

A Markov chain is oompletely described if the initial 

probability function (pl^^ is the probability that 

the chain is in state i at time zero) , the one-step 

transition probabilities and the state space of the 

chain are given. A finite Markov chain is ergodic if 

there exist probabilities n. such that - ir 

3 n-*-« 3 

for all i and j. is the probability of going from 
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state i to state j in exactly n steps. All states of the 
Markov claain belong to tie sane ergodio class, if there 
exists an n such that > 0 for all j and k. ¥6 have a 

a theorsn governing the limiting behavioux* of an ergodic 
Markov chain with all its states belonging to the same 
class by which we have 


S 

(i) > 0 and 5! 

k=1 



S 


O 


(1.3.1) 

(1.3.2) 


¥e shall use these results in arriving at a Monte Carlo 
scneme for the estimation of theraodynamic properties of 
a classical system. 


Squilibriu-i statistical aschanics is concerned with 
averages of the type 

' ' << 

over all over all 

states -BE. states -SS . 

<5^ = I ® /I e (1.3.3) 

i i 

where is the energy of state i, is the value of the 
variable F in state i and B = * This sum is reduced 

to an integration over configuration space for classical 
systems. ¥e would however be concerned mostly with 
properties of the system which depend only on space cooordi- 

nates . Hence the above sum would be reduced to an inte- 

*■ 

gration of the form. 
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<F> = jT{r) dr // dr (Ii3»^) 

■where U(r) is the potential energy of interaction among the 
particles and r designates all the space variables of the 
system. Thus evaluation of average properties of the 
system reduces to an evaluation of aulti-dimensional 
integrals. The exponential factor in the integrand implies 
that a significant part of the integral is concentrated in 
a very small region. Hence estimation of <F> appears 
at the outset to be Just a straight -f orward estimation of 
the multi-dimensional integral by an 'importance sampling' 
method in -which the random variables are drawn from a 
distribution given by dr. This means that 

we should know the denominator to get the distribution 
function. Ihiowing the denominator is equivalent to solving 
the problem itself since it is proportional to the partition 
function of the system. We set out to estimate the thermo- 
dynamic properties hy Monte Carlo method since detenaina- 
tion of the partition function by other means is not pos- 
sible. Now it appears that we have ended up in a vicious 

circle. A way out of this was devised by Metropolis, 

• 35 ■ - 

Rosenbluth, Eosenbluth, Teller and Teller. They devised 
a Markov chain such that the ratio of the limiting prot^- 
bilities of any two states of the chain is equal to the^^^^^^^ ^ 
ratio of the equilibrium probabilities of the corresponding 
two states of the system; 
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( 1 . 3 . 5 ) 


Here, it. and. ir. are the limiting probabilities of the 
states i and j of the Markov chain. Hence if the Markov 
chain is propagated for a sufficient length of time, the 
states of the Markov chain would be distributed with 
probabilities given by the equilibrium probabilities of 
the given system. Thermodynamic averages of the system 
can now be estimated as the average over all states of 
the Markov chain, taken over a large number of states, 
after the Markov chain has attained limiting behaviour. 

In this method, since only the ratios of the limiting 
probabilities of any two states need be known, the problem 
of evaluating the denominator in eq. (l.3.^) is overcome. 

35 

The Markov chain of Metropolis et al. for a system 
under study is defined by initial probabilities taken to be 
the thermodynamic probabilities for the states of the 
system in equilibrium at T=0 or T= «> (completely ordered 
or' completely disordered system) . The one step probabili- 
ties for the chain to go from a state 'r* to state 's’ is 
unity if the energies of the corresponding states of the 
system are such that E > E otheiwise it is given by 
exp [(E r® )sl . The state spsMse of the Markov <dialn 
cbn^jletely defined and the process of going from one state 
to other is such that all the states of the state space 
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are accessible from any state. These conditions ensure 
that corresponding to any two states 'r' and 's' of the 
canonical ensemble, we generate v and v„ systems, 

X* s 

respectively, such that 




)S 


{ 1 . 3 . 6 } 


This is because the probability with which we pick a state 

'r< to move to state 's' is equal to the probability with 

which v/e pick the state 's' to move to state 'r' , i,e., 

P = F , Now , if S > S , the number of systems moving 
rs sr ’ r — s 

from state 'r' to state 's' is v P and the number of 

r rs 

systems moving from state 's' to state 'r' is v P 

s rs 

exp [-(B -B )6]. The net number of systems moving from 

X* S ' 

state 's' to state 'r' is, 


N = V P '.exp [-(B -B )6] - V ? . (1.3.7) 

If < exp( -E^6)/exp{-E^ g) on an avereige more 

systems move from 's' to 'r' . This coupled with the 
condition that all states should be attainable by the 
process enables the relation (1.3.5) to be satisfied. 

» 

The one step probabilities .assumed by Metropolis 
35 ' 

et al. to define the Markov chain is not unique. ¥ood 
and Parker^ recognized that the Markov chain we are^^^ 
dealing with is an ergodic Markov chain with all its states 
belonging to the same class and used the relations (1.3.3) 
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and { 1 . 3 * 2 ) to arrive at the conditions the one step 
probabilities of such a Markov chain should satisfy* By 
inspection » we see the conditions are normalization and 
microscopic reversibility, i.e. 


I = 1 for all j = 1,2,.. .3 and (I. 3 . 8 ) 

k=1 




= ^k^kj- 


(I -3.9) 


In our case the limiting probability of the Markov chain 

should be proportional to the thermodynamic probability 

at equilibrium, ir. = c exp (-E.g). ¥ood and Parker used 

D 1 

the following one step probabilities to obtain Monte Carlo 
equation of state of interacting molecules 




and p . . 


exp F-{S^.-Sj,)e] 


for E . > E, 
3 Ic 


(I. 3 . TO) 


= 1-1 
k;^j 




ak 


Computation methods for Monte Carlo estimation of 
thermodynamic properties is straight forwar^. The system 
'is created on the computer* with the given initial 
probabilities. The location (inside the system) at which 
a change is to take place is selected by a random process. 

The nature of the change at the location may be probabilistic 

*The present study wias carried out oh an 134 7044/l4ol 
coE^uter'system of this Institute. 
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or deterministic, depending on the problem. In the jformer 
case the nature of the change is decided with the help of 
a random number drawn from an appropriate distribution. 
Probability for the change is computed from the one step 
probability functions. If the probability is unity, the 
change is carried out. If the probability is less than 
tmity, a random number between 0 and 1 is generated and 
the change is effected if the random number is less than 
the probability, if not, the original state is retained. 

After repeating this process sufficient number of times, 
the effect of initial configuration vanishes. From this 
point onwaiKis, in the calculation, exact thermodjnaami c 
variable corresponding to the state of the system is 
calculated at each step. The ensemble average is taken 
to be the average of the thermodynamic properties over 
all the Monte Carlo cycles perfomed after we st anted 
oon^juting the thermodynamic variable . Unsuccessful tries 
should also be included in the averaging process as not 
to include them is equivalent to removing from the 
ensemble those systems which tried to move but were 
forbidden. 

In all Monte Carlo calculations, it is isiportant to 
know, (i) the statistical error^*^’^® involved in the method., 

and (ii) the number of Monte Carlo steps required for the^ ^ 

' 37 

effect of initial oonf iguration to vanish, thereby bringing 
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the system to tiiermodsnaami c equilibrium at tbe specix'ied 
values of the external variables (temperature, magnetic 
field etc,). To estimate these, let us start with the 
relations defining the relaxation fvmctions 

<^jg(t) = {<A(o) B{t)> - <A><B>>/ {<AB>-<A><B> 

(-fc) = 

{<A(0)>^, - <A(»)>^,>, (1.3.12) 


A and 3 are the variables of inters s *• hASB 
librium relaxation function and gives the dynamics of 
fluctuations around equilibrium. 4*^^ (t) is the non- 
equilibrium relaxation function and gives the dynamics 
of change in A as the external variable e is changed from 
e’ for t < 0 to e' + Ae at t=0. The averages in the 
expression are calculated over the states for- which the 

external variable is e' , Using these functions two 

AG * 

relaxation times, defined by the 

relations 


"5A6B = ^ VSB (*) 

and (t) at. (1.3.1'*) 

The number of Monte Carlo steps, n, that should be excluded 
to get rid of the effect of initial configuration, after 
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21 ^ i© , 

the chajige Ae is effected, must satisfy wnere 

N is the size of the system. The expected statistical 
error iSA of an estimate A, after m Monte Carlo steps is 


(6A)^ 'V ^ [ <A^> - <A>^3 (1 + 2T 


V m 

)t - »T 


6A5A^* N 


5A5A * 


The quant it i e s x 


Ae 

6A 


T . can be estimated a priori 
AAA A 


from the theory of stochastic model systems . 


In the present study, instead of making these esti- 
mates, "we have relied on comparison with exact results for 

testing the accuracy of the method. Such comparisons have 

38 39 

been made by earlier workers ’ in other statistical 
mechanical problems and Honte Carlo method has been shown 
to be quite accurate. The main sources of error in the 
Monte Carlo method are due to finite sample size, finite 
size of the system and systematic errors. Finite sample 
size leads to (a) a statistical error, proportional to 
thermodynamic fluctuations, which is always present even 
if the san^le is drawn from a true ensemble with perfect 
randomness, and (b) an error present due to correlation 
among states closely succeeding one another. The error 
from (b) has the effect of marginally reducing the efficiency 
of the Monte Carlo method. Errors due to finite size of 
the system can be eliminated using theories which provide 
in foniKit ion about the finite size effects. In a study of 
phase transitions finite size has the effect of smoothening 
out the variations of appropriate the rmodsmamic quantities. 
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which othen^ise should show sharp breaks at th© transit ion 
points. Systematic errors can be eliminated by using good 
random variables in the calculation and by proper debugging 
of the programs . 

Monts Carlo method has been applied in statistical 

mechanics to study a variety of* problems such as order- 

^0 

disorder transformation in b.c.c. systems, phase tfefl3.ai- 

tion of* three and two-dimensional Ising system with a variety 

-U-3 

of interactions, equations of states of molecules 

35 35 44 

interacting with different potentials ' ’ and so on. 

35 

The method of Metropolis et al. is applicable only 
to obtain canonical ensemble averages of classical systems. 

A study of phase transitions by canonical ensemble method 
requires a fairly large system. This is not a problem in 
one-dimension; However in 2 and 3-<ii®®iisions the system 
size required in some cases may become unmanageably large. 

¥e can overcome this difficulty if we can obtain thermo- 
dynamic averages of a grand canonical ensemble by Monte 
Carlo method. The system size required to study a phase 

transition by dealing with a grand canonical ©nsemble is 

45 

quite small. Pi lino v and Norman have given a method by 
which we can generate a Markov chain with state space 
consisting of the states of a grand canonical ©nsemble . 

To deal with quantum statistical mechanics by the 
Mont© Carlo method, we have to make modif icat ions of the 
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method of* Metropolis et al. The problem arises since 

in quantum statistical mechanics the partition function 

cannot be represented by a conf i^rat ional integral and 

we need to know the eigen values of the system, which is 

46 4 ? 

often diffidult to obtain, Handscomb and FoSdick have 
given methods for estimating thermodynamic properties of 
the system by a Monts Carlo calculation. Handscomb»s 
method involves a partitioning of the Hamiltonian and 
expressing the thermodynamic averages as a sum of traces 
taken over different parts of the Hamiltonian and estima- 
ting this sum by a. method similar to Metropolic method, 

provided the traces appearing in the sum are easily 

4 ? 

calculated,. Fosdiok writes the partition fimction as a 
Wiener integral and represents this Wiener integral by 

f 

an n-fold Riemann integral, the evaluation of the latter 
is done by standard Monte Carlo methods . 

Monte Carlo method has been applied for the extrac- 

48 

tion of critical exponents of a system. In this method, 
the transformat ion of interaction parameters, defining 
the model, imder a change of length scale is obtained by 
Monte Carlo method . The avereige lifetimes of a single spin 
in up and down state and the average lifetimes of a block 
of spins in up and down states of the block, are estimated 
by Monte Cairlo method. These give the transformation of 
parameters under a change of length scale, which in turn 
is related to the critical exponents of the system. 
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1.4 Programs and Testiiag 

In the present study, we have simulated an Ising 
■ch^iin. 'With one thousand spins. The up and down states of 
a spin are designated by the numbers one and. two. To begin 
with, the interaction parameters, the temperature and a 
variable which decides the initial state of the chain are 
read in. If the initial state to be generated corresponds 
to a completely ordered chain, number 1 is assigned to all 
the sites. If the initial state to be generated happens 
to be a disordered state, we pick the first site, call a 
randoni number between 0 and 1 and assign 1 to the site if 
the random number is less than 0.5 and assign 2 otherwise, 
on repeating this with all the sites we get the chain in 
a disordered state. The initial energy of the chain is 
con^juted in two parts viz., (i) the contribution from 
short-range interactions and (ii) the contribution from 
the long-range part. The contribution from (i) is calcu- 
lated by starting from the first site and evaluating its 
interaction energy by finding out the state of the spins 
at the first and second -neighbour positions. This is 
repeated for all the thousand sites and the total ^ort- 
range interaction energy is obtained by summing up all 
these contributions. A cyclic boundary condition is used 
in all the calculations. The long-range contribution to 
initial energy is given by '^^r ^ where 


is the long-range 
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interaction strength and M is the magnetization of the 
ohaTn in the initial state. Magnetization is given by the 
number of '1 sites' minus the number of '2 sites'. In all 
our calculations we have used equivalent -neighbour inter- 
action for the long-range part. The initial energy is the 
sum of contributions from (i) and (ii). 

The chain, which is now at T=0 or T=«, depending on 
whether its state is ordered or disordered, has to be 
brought into equilibrium at the temperature at which we 
intend to calculate the thermodynamic properties. This 
is achieved by flipping the spins in the followong way, 

¥e select a site at random with the help of a random 
integer between O and 1000, ¥e calculate the change in 
energy of the chain, AE, when the spin at this site is 
flipped. &E is calculated as the sum of ^^sr' change 

in short-range interaction energy of a segment'’ of the chain 
before and after flipping the spin (the segment of the 
chain includes all the spins whose short-range interaction 
siicsgies are affected by the spin flip) and change 

in long-range interaction energy caused by the charge in 
magnetization accon^janying the spin flip. The spin is 
flipped if AE is negative or zero. If AE is positive, 
a random number between 0 and 1 is generated and con^jared 
with P given by exp(-AE/T)« If the random number is 
smaller than P the spin is flipped otherwise the original 
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state is retained. The total energy of the chain is 
changed by AE if the spin has been flipped otherwise the 
original value is left unchanged. The whole process is 
Treated starting from the selection of a random site. 
Depending upon the temperature, 3000 to 6000 O tries are 
required to bring the chain to equilibrium. In equili- 
brium, the intemal energy starts fluctuating around a 
mean and the test of having reached equilibrium lies in 
looking for fluctuations in the internal energy which is 
printed out at the end of every hxmdred or hundred and 
fifty spin flips. 

Ones having ensured that the system has reached 
equilibrium, we can start con^juting the thermodynamic 
properties. The thermodynamic properties computed are 
the internal energy, square of the internal energy, 
sublattice magnetizations and correlation function 
describing correlation of spins on different sites. To get 
thase , the respective sums (formed by adding the actual 
value of the quantities at the end of each Monte Carlo step 
including unsuccessful tries) are calculated. The avereiges 
are given by dividing these sums by the number of tries 
made after first coq^jutation of these quantities. The 
sublattice magnetizations are estimated asstzming the 
existence of two sublattices viz. all the odd numbered 
sites fonaing one sublattice and all the even numbered 



sites fonaing the other. The correlation fiinction of the 
chain is initially calculated as. 




y. v(i) v(i) 

j 3 0+k’ 


k = 1, 2, 


...20 


(1.4.1) 


where (i) refers to the initial calculation. Subsequent 
calculations of g(k) involved making appropriate changes 
whenever a spin was flipped. Specific heat is given by 
fluctuations in internal energy by the relation 


C {<E^> - <S>^> (1.4.2) 

2 2 

Since we have cornputed <E > as well as <S> , specific 
heat estimate can be maide from our calculations. Uniform 
susceptibility is given by 

X a f fj g(r,r') dr dr» (1.4.3) 

where g(r,'r') is the correlation function for magnetization, 
which we have estimated for our discrete chain and hence 
we can calculate susceptibility by our simulation. In 
fact, we have estimated staggered susceptibilities, 
assuming different orderings of the spin, by forming 
appropriate partial sums of the correlation fmction. 

In all our calculations the number of tries requlaTed 
to reaoh equilibrium at atny given ten^erature is read in at 
the beginning of the program. This is estimated, over all 
tea5>eratu37es of interest, from trial runs on a puirely anti- 
ferromagnetic chain with a disordered initial configuration. 
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An examination of the internal energy that is printed out 
at intervals of every htandred spin flips shows that the 
number of tries excluded from the averaging process is 
generally sufficient to remove the effect of initial 
configuration. Averages of thermodynamic properties are 
taiken over sixty to seventy thousand tries. The time 
taken per Monte Carlo cycle was between 3 to 6 milliseconds, 
depending on the temperature. At the end of the Monte Carlo 
calculation, the energy of the chain is calculated ab initio 
and coiiT>ared with the energy of the chain at the end of 
the Monte Carlo cycles. These two should agree to within 
the round off error and this forms one of the checks on 
the correctness of the calculation. Other checks of the 
program included tracking of all the variables involved 
in the calculation by printing them out as and when they ap- 
peared in the calculations. 

Accuracy of the calculations is checked by performing 
the calculation on a nearest -neighbour antiferromagnetic 
Xsing chain over the teo^erature range of interest. The 
exact thermodynamic properties can be easily evaluated 
from the partition function which can be determined easily. 
Fig. 1.2 gives the Monte Carlo estimation of specific heats 
and internal energy as well as the exact specific heats 
and internal energy, over eighteen tea^eratures in the 
teagjerature region of interest. The agreement between^^^^^^^^^^ ^ 
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Monte Caxlo and exact values is good in the higher temp- 
erature region. But the peaks in the specific heat curves 
from the Monte Carlo method and the exact calculation' do 
not appear at the same teoperature . This could be due 
to the presence of large domains, in the Ising chain 
simulated in our calculations, at low temperatures. The 
specific heat estimates in the temperature region 0 to 
T “ 2j/3k is not reliable because the number of breaks in 
the chain ordering are far in excess of the equilibrium 
number of breaks expected at that temperature. This cannot 
be got rid of as the number of Monte Carlo cycles required 
to remove a single break is rather large. This problem 
may be overcome by using different spin-flip mechanisms, 
which would make the chain ergodic at these temperatures 
also.' However, from the work of Nagle and Theumann and 
H/rfye we knovj that this temperature region is not of interest. 
The sublattice magnetization computed in our calculations 
do not agree with the magnetization from the exact calcula- 
tions for the same reason, viz. existence of domain structxare 
Pig. 1»3 gives the correlation function of the antiferro- 
magnetic chain at three different temperatures. 

¥e have performed all our calculations by assuming 
an initially disoiPdered state. In the case of nearest- 
neighbour antiferromagnetic chain we see that the results 
in the high teoperfiture region do not depend ipon the ihit i a 1 
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state. Since our interest lies also in ini ei^ne ting 
polytypism where the crystals are grown at high teapera- 
tures, ws decided to carry out all our calculations using 
an initially disoixiered state* 


1*5 Results and Discussion 

The values o£ interaction parameters for which the 
calculations have been carried out are listed in Table I.l. 


Table I 


Values of 
which 

interaction 
the present 

strengths and 
calculations w 

magnetic fields for 
ere carried out. 

Sr .No. 

H 

J 

K 

A 

1 

C.O 

12.5 

0.0 

30.0 

2 

2.0 

12.5 

0.0 

30.0 

3 

0.0 

6.0 

6.0 

30.0 

k 

1.0 

6,0 

6.0 

30.0 

5 

0.0 

6.0 

12,0 

30.0 

6 

2,0 

6.0 

12.0 

30.0 

7 

0.0 

12^5 

12.5 

30.0 

8 

4.0 

12.5 

12.5 

30.0 

9 

0.0 

15.0 

9.0 

30.0 

10 

3.0 

15.0 

9.0 

30.0 

11 

0.0 

18.0 

7.0 

30.0 

12 

3.0 

18.0 

7.0 

30. 0 

13 

G.O 

15.0 

3.0 

30.0 

14 

2.0 

15.0 

3.0 

30.0 


H = Magnetic field, J = Nearest neighbour interaction 
strength, E s= Next nearest neighbour interaction strength 
and A/N = A/1 000 s= Equivalent neighbour interaction strength. 



37 


Calculations were performed botli in zero field and in ' > 

the presence of a field. The first two calculations 
corre spend to Nagle’s model and the rest are based on 
Theumann and H/iye»s model. 

The values assumed by the interaction parameters ^ 
in our study of Nagle’s model in zero field correspond to 
those values for which a first order transition is predic- 
ted by Nagle. Our calculations do not show the oocxirrence 
of a first order phase transition in specific heat vs 
ten^jerature plots (Fig. 1.4) and the phase is essentially 
antiferromagnetic at all ten^jeratures. However, in the 
presence of a field, we observe a first older phase transi- 
tion from a low -temperature long-range ferromagnetic 
ordering to a high-temperature short-range ant if erromagne tic 
ordering, ¥e do not see the predicted first order phase 
transition in zero field apparently due to the fact that 
the system is in a metastable state below the critical 
ten^xerature . The low-ten 5 >erature phase in zero field that 
we observe is ant if erromagnetic while the ground state 
internal energy is lower for a ferromagnetic phase. Even 
in our calculations in the presence of a field, the low- 
teix5>erature phase is antifenromagnetic while once again 
the ferromagnetic phase is more stable at T=0 . In f act , 
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at one of the tesaperatures, in the presence of a field, 
we observe the system going from a metastable state to 
a stable stats. Fig. 1.5 gives the internal ®aergy as a 
function of the number of Monte Carlo cycles for this 
case. The antiferromagnetic state satisfies all the 
criteria that should be satisfied by a metastable state: 

(i) presence of a single phase (in this case, the antifer- 
romagnetic phase) through out the system, (ii) very long 

life-time of the state as evident from the figure, and 

* 

(iii) very little likelihood of returning to the metastable 
state once the system has reached a stable state. The 
metastable stats decays into a stable state ovdr rela- 
tively few Monte Carlo cycles, as expected from the 
kinetics of this decay. 

IfQ 

Lebcxwitz and Penrose showed that for systems 
interactdng: with a short -ran^e potential and an inf inite- 
rang© Kac potential, there exists at least one region 
in configuration space such that the escape rate for the 
system from that region of configuration space is arbitra- 
rily small as also the equilibrium probability for the 
system to be in that region of configui*ation space* 

They further shewed that if the system enters a metastabie 
state, then phase transitions arising from the infinite- 
range Kao potential do not manifest themselves* However , 
in higher dimensions, the phase transitions arising due 
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to sliort - 3 ? aixge potentials , are observed independent of 
whetlier the system has entered a metastable state or not* 
Our system being one**dimensional, phase transitions can 
occur only due to infinite -range Kac potential and in the 
present study > the 33^3 tem having Entered a metast able state 
fails to show the phase transition. Binder and Stoll^^ 
have also observed such metastable states in their study 
of square lattices. They find that the time taken for the 
domains in a ferromagnet to reorient themselves with the 
field, when the direction of the field is reversed, is 
much longer than defined by (l.^,lU), Experimentally 

also cases like, magnetization of a domain (in a ferro— 
magnet) remaining opposed to the direction of an applied 
field, are well-known. Our calculations on Theumaim 
H/$ye*s model also showed such metastable states in some 
regions. We can possibly avoid such metastable states in 
simulations by starting with a different initial state*. 

Our calculations on the Theuraann-H^ye model are 
aimed at understanding the way this Ising chain behaves 
as a fimction of temperature in various regions. It is 
expected to show several phase transitions as a function 
"fe ©mperat ure ( th e ntamber of trans it ions dep ending upon 
the regioa of the phase diagram shown in Pig. 1 . 1 ) accom- 
panied by changes in the ordering of the spins. 

^ \ ' • ' . ■ • 

*W© have not carried out this calculation since our interest 
lies mainly in the Theumann-^^^ model. 
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In region I, the system is in a low-tes 5 >eI*at^l^o 
ferromagnetic state. With increasing teisperature , the 
chains show a first order transition with a large change 
in internal energy and a break in specific heat cui^e 
(Fig. 1.6), The- correlation f motion for various teBq>era- 
tures is shown in Fig. 1.7. The high -temperature phase 
shows a short-range ordering of spins described by a unit 
y cell with 4 lattice sites having the spin configuration ++>+ . 

Jn the presence of a magnetic field, we observe a metsstable 
state at low -temperatures in which the magnetization of the 
chain and the applied field are in opposite directions. At 
higher temperatures, the stable state of the chain is 
restored and this is obsearved to xmdergo a first order 
phase transition. Once again this first older phase 
transition is accompanied by breaks in internal energy 
and specific heat curves (Fig. 1.6). 

pur calculations in other regions do not show such 
sharp phase transitions as in region I. Instead, we observe 
small peaks in specific heat curves. The intemal energy 
curves at these points are no longer smooth. The short- 
range ordering of the chain in each of these regions does 
not show any significant change with ten^terature althoiagh 
the low t,ea 5 )erature oidering of the chain is quite different fo 
different regions. For the values of interaction parav^eter 
falling in region II, the specific heat versus temperattjre 
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curve (Fig. 1.8) shows two peaks indicating the existence 
of two phase transitions of second or higher order. In the 
presence of a field, one of the peaks veinishes indicating 
that the magnetic field has suppressed one of the phase 
transitions. In region III, the specific heat versus 
tenj>ea?ature cuive shows as many as four peaks showing that 
there are four pheise transitions for the system with that 
set of interaction parameters. These peaks disappear in 
the presence of a field showing that in this case also the 
phase transitions are suppressed by a magnetic field. The 
specific heat curve in region IV does not have any clear 
peaks though there appear to be some broad huo^s . In the 
presence of a field, these humps disappear except the one in 
the higher temperature side. This becomes more pronounced 
giving indication of a phase transition in the presence of 
a field. Similarly, in region V, the specific heat oxarve 
shows only broad hun^s which disappear in the presence of 
a field. In region VI, the specific heat curve shows two 
prominent peaks showing the occurrence of two phase transi- 
tions. In a msgnetic field the peaks shift to the higher 
teH^^erature side in such a way that the two peaks come close 
together. 

At low teH:^ei*atures, the chain exhibits different 
kinds of ordering in . the different regions. Tt® ordering 
of spins in region II at low temperatures is described by a 
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unit cell of four lattice sites with spin cooof Iguration 
in the unit cell being ++++ . In region III, the chain 

has some segments over which the ordering can be described 
by the \mit cell with four lattice sites and over the 
remaining segments, the ordering is given by a unit cell 
of 3 lattice sites per cell, in the spin configuration '!•++ 
or <.+4. » The chain in region IV is having a spin ordering 
given by a unit cell with 3 lattice sites with the confi- 
guration in the cell being +4.4.* The ordering in region V 
is essentially the same as in region IV except for small 
stretches of an ti ferromagnetic oi*dering in the chain. 

The ground state of region VI is described by unit cells 
with 3 lattice sites per cell (spin configuration in the 
cell can be either +++ or 4-++ ) extending over short 

stretches in an otherwise antiferromagnetic chain. 

Correlation fxmotions, up to eighteenth neighbour, 
of the chain for regions II, III and IV are shown in 
Figs. 1.9 to 1 . 11 , for scHse temperatures. ¥e see that, 
giyen the region, the short-range oiviering remains essen- 
tially the same over all temperatures and no distinct 
changes are brought about by phase transitions. Magnetic 
field does not change the orderings significantly, thoiigh 
the unit cells with 3 lattice sites tend to have their 
magnetization oriented along the field. 
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Sub lattice magnetizations computed assuming diflTe rent 
unit cells over all regions, are not of mu ch help in finding 
out the orderings of a chain at different temperatures. 

This is due to cancellation of these magnetizations caused 
by the domain structure of the chain. Hence, direct computa- 
tion of sublattice magnetization is dLncapable of reflecting 
the actual ordering of the chain. This problem can be 
overcome by computing sublattice magnetizations from the 
correlation functions. This is possible since given the 
spin orientation at a site, the correlation function., gives 
the probability that the spins at 1st, 2nd, ... neighbours 
also have the same orientation and the domain stnicture 
does not mask the correlation function from reflecting 
the ordering in the system. Thus, for example, by su mm ing 
up the value the correlation function assumes at 33C’d, 6th, 
9th, ... neighbours, we get sublattice magnetization 
corresponding to a 3 sites per tmit cell structure, after 
proper normalizations. ¥e have obtained sublattice magne- 
tization corresponding to different spin orderings . 


In Fig. 1.12 we have shown different sublattice 
magnetizations as a function of ten^erature for region II. 

We do not see ainy significant changes in oidering as a 
function of temperature. If there wore to be different 
Spin ordering s , at different ten^eratures, at phase transi- 
tions , we should have observed some s^tjl^'^jgl^.^jj^petization 
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abmaptly decreasing while some other sublattice magnetization 
increase in magnitude. The reason why we fail to see such 
changes is that the short-range order continues to be the 
same over all ten 5 >eratures. And, even if the long-range 
order is different, the oontribrtion from the long-range 
order terms is very small and hence is not seen in the 
sublattics magnetization curves. 

The long-range order of the system can be obtained 
from the asyn5>totic correlation functions. ¥e con 5 >uted 
correlation functions for region II (this region was chosen 
because there are only two phase transitions in the region 
as seen in Fig, 1,8 and hence is the simplest to deal with) 
starting from 23 rd neighbour up to 4oth neighbour, over all 
tec^eratures. These are shown (Fig. 1.13) for different 
ten^eratures, at which different phases are expected to 
exist. For comparison asjnmptotic correlation functions 
at these ten5)eratures are given for ths nearest -neighbour 
antiferromagnetic chain (Fig. I.l4). ¥e see that at low 
teii^eratures, the long-range order is described by a unit 
cell with 4 lattice sites ( ++++ ) • ¥ith increase in ten^)- 
erature, the system goes over to a phase where the long- 
range order is now described by a unit cell with 3 lattice 
sites ( ). After the second phase transition, the system 

tends to have a ferrcraagnetic alignment, though it is not 
seen definitely. These results sre in keeping with the 
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long-range ferromagnetic interaction term present in the 
Hamiltonian. ¥e see the system racx'iig over from a zero — 
magnetization phase to a magnetically ordered phase in 
two steps, with increasing ten 5 >erature , as is expected from th 
loaager range of the ferromagnetic interaction. The 
asyn^jtotic correlation functions also reflect the domain 
sizes. The distance over which the anplitude of tlw 
asyrspt otic correlation function remains roughly constant 
is a measure of the size of domains in the system. 

The present calculations establish that the Theumann- 
Hjiiye model exhibits phase transitions as a function of 
tenperature , as predicted. These phase transitions are 
accompanied by a change in long-range ordering of the 
spins although short-range ordering remains essentially 
the same over all temperatures. The number of phase 
transitions in a given region as a function of temperature 
appears to be the same as the ntunber of transitions at 0 Z. 

(in that region) as a function of interaction strengths, 
Theumann and H;^ye predicted first order transition in all 
the regions, at OK. Our calculations, however, do not 
show first order transition, as a function of tenperature, 
in any of the region s except region I. 
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1.6 Summary 

Spin-half Ising chains with short-range antiferro- 
aagnetic and infinite -range ferromagnetic interactions are 
known to exhibit several phase transitions, Nagle observed 
a first order transition in the presence of a field as well 
as in zero field as a function of te:i 5 )e nature when the 
short-range interaction was restricted to first neighbotir 
and the infinite range interaction was talcen to be an 
equivalent neighbour interaction. Theumann and H/ye 
extended the short-range interaction up to second neighbour 
and used infinite-range Uac potential for the long-range 
part. They obtained a phase diagram with interaction 
strengths as variables for zero temperatui^ and found that 
in some regions of this phase diagram, the chain showed 
as pany as four first order phase transitions . 

In the present study, we have carried out Monte Carlo 
calculations on the Ising chains of Nagle and Theumann and 
Tifiye . From our calculations we have obtained estimates of 
specific heats, internal energies and correlation functions, 
as a fxmetion of temperature for various interaction streng- 
ths for both the chains, m the case of Nagle 's Ising chsdLn, 
wo have observed the predicted first oinier phase transition 
in the presence of a field, ¥e did not, however, see the 
predicted first order transition in zero field and we have 
attributed this to the formation of a metastable state at 
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low ■feeB5)eratures. In the case ot ThBuma.im-lifiyQ Ising chain, 
in one of the regions of their phase diagram, we have 
obsei*ved a first order phase transition as a function of 
tesi^erature , both in the presence and absence of a field. 
This was accoB^janied by a change in the short-range 
ordering of spins. In the other regions of the phase 
diagram, wo did not see such first order phase transitions, 
as a function of ten^erature, although the specific heat 
vs temperature curves showed distinct anomalies indicating 
occurrence of second or higher order transitions. The low - 
tejEperature spin ordering of the chain is found to be 
different in different regions of the phase diagram . The 
specific heat anomalies were not accompanied by any percep- 
tible change in the ordering of the spins, as seen from 
the microstates. Sublattice magnetizations, (computed from 
the correlation fianctions) also do not show any characteris- 
tic changes as a function of temperature. Eowever, asymp- 
totic correlation ftinctions of region II showed that the 
long-range ordering of the spins changed when the chain 
landerwent a phase transition. The different long-range 
oixiering of spins in different phases was explained on the 
basis of longer range of ferromagnetic interactions. In 
addition to providing support to the models of Nagle and 
Theumann and ye, the present Monte Carlo study provides 
details of the phase transitions and spin orderings as a 
function oft emperat ure in various regions of their piiase 
diagrams. 
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Conq>lete characterization of the phases in each 
region, as a ftmction of temperature, is plausible if a 
free energy expansion in termis of various order parameters 
can be worked out. One possible method of arriving at 
such free energy expansions seems to be the fxmctional 
integral method. ¥e propose to carry out such calcula- 
tions in the future to understand the changes in spin 
ordering accompanying these phase transitions. 
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APPENDUC 

3AND0M NUM3S3. 3E1-1SHATI0N AlTD TESTING 

In this, we otitlins tbs netliod of randoc nartbsr 
gsnaration and the tests performed to ascsrtaiii that thsy 
belong to a uniforrn distribution. .?,andoE numbers were 
generated by th® following algcritha. An integer, ini- 
tially chosan to be 377773 {i^ octal) was multiplied by 
(2 +3). If the resulting integer is less than 2^ -3, 

it is taksii as the first random number. If it is greater 
than 2'^'^-1 (this number is specific to IHi 7044 on which 
all our calculations were carried out) , the computer 

truncates the digits of t3"e number from the left till the 

35 

resulting number is less than 2-1. The resulting number 

is tahsn as tbs random number. For the next cycle, the 

last random number that is generated is multiplied by 

(2 V+3) and the whole process is repeated. Handom numbers 

between 0 and 1 are obtained by di'viding the above random 
35 

numbers by 2 . Only two operations are required to 

generate a random number by this method. 

The first thirty thousand random numbers generated 
by this method were subjected to the following tests to 
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ascertain that they belonged to a uniform distribution ^ 

(i) Mom ant 5 test, where up to fourth moment of these 
numbers were calculated and coii^ared with the expected 
moments* (ii) Frequency test, where the actual frequency 
f^ of the random numbers is calculated and con^ared. with 
expected values* f^ denotes the number of random numbers, 
from among the thirty thousand generated , that belong to 
the ith sub-inteival defined by the end points (r-l)/10 
and dL/lO* (iii) The correlation tests in which the fre- 
quencies are coinputsd* Here i denotes the ::,th sub-j. 

interval to which a random number belongs and j the jth 
sub-interval to which the next random number belongs • The 

superscript (1) implies that v/e are testing for correla- 

2 

tion between numbers generated in succession# X -test 

is not reliable in these cases as the division of the 

2 

interval (0-l) is arbitrary* However, the X —test for 
both correlation and distribution show that the random 
numbers generated belong to a good uniform distribution* 
Results of the tests (i), (ii) and (iii) are given in 
Tables lA, IIA and IIIA respectively . 
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Tabl* IDL 
Moments Test 

(Sxpeeted values are calculated for uniform distribution) 


Order of 
the moment 


2 3 



Actual 0.5003 0.33^3 0.2513 0.2014 

value 


Expected 0.5000 0.3333 0.2500 0.2000 

value 


Table IIA 
Frequency Test 


Interval 

0-.1 

r -A I 

.1-.2 

.2-. 3 

.3-4 

. 4 -. 5 

. 5-. 6 

.6-. 7 . 7-.8 

V ■ r 

. 8 ^ .9 

.9-1 *0 

Actual 

frequency 

3106 

2935 

2977 

3021 

2937 

3045 

2918 

2875 

3189 

2996 

Expected 

frequency 

3000 

3000 

3000 

3000 

3000 

3000 

3000 

3000 

3000 

3000 
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Table IIIA 

Frequencies for correlation of successive random numbers 
(Expected frequency: 300 for all cases) 



- - 




_ " i-fc. A. A' 


-- - _ A,_ *. 




i-*- / 

r ^ * 

0-^1 

. 1-.2 

. 2-. 3 

« 

1 

* 

. 4-. 5 

.5 -. 6 

. 6-. 7 

.7 -. 8 

.8-. 9 

.9-1 

0-.1 

356 

329 

297 

322 

275 ' 

313 

317 

257 

337 

303 

. 1 - ,2 

267 

290 

314 

281 

295 

310 

321 

280 

300 

277 

. 2-. 3 

303 

288 

295 

318 

311 

3 o 4 

249 

287 

318 

3 o 4 

. 3 -.k 

316 

310 

306 

274 

301 

309 

280 

299 

310 

316 

A -,5 

298 

268 

284 

306 

306 

292 

300 

256 

325 

302 

. 5 -.^ 

337 

295 

288 

299 

311 

297 

296 

307 

306 

309 

. 6-. 7 

319 

283 

286 

278 

271 

316 

288 

273 

295 

309 

00 

• 

1 

• 

292 

284 

273 

287 

262 

288 

288 

3 o 4 

337 

260 

.G -.9 

327 

287 

332 

321 

315 

320 

305 

336 

337 

309 

• 9 - 1 .0 

291 

301 

302 

335 

289 

296 

274 

276 

325 

307 

: 

_ lA - -k. 

1 - r 
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CHApraR II 


MCaSTS QARIiO SIMUIATIOH CP POLYTYPES 


II . 1 pat^^dact ion 

Close packed structures ir. solids can be thought 
of as obtained by packing atonic layers, one above another 
Two connonly knov^n close packed structures are the hexa- 
gonal close-packed and the cubic close -packed structures. 
Hexagonal close -packed structures are obtained by packing 
layers of atons such that atoms in alternate layers come, 
exactly, one above another. If we can represent the 
position of atoms in a layer by the alphabets A and B, in 
hexagonal doss packing the layers are stacked in the 
sequence . . .ABABAB . . . . . In cubic close packing, the atoms 
of every first and fourth layer are superposed on one 
another. In the above notation, the layer sequence would 
be .. .ABGABC... This is because, after the first two 
layers are packed one above another, ws eire left with 
only two options - one leading to a cubic close packing 
and the other to a hexagonal close packing. A given 
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atomic layer is said to be in a hexagonal close-pecked, 
configuration (h) if atoms in its neighbouring layers 
are exactly superposed on cm *sii>thor(6.g., layer 'B* in 
the sequence ...ABk,..). Otherwise, the layer is said 
to be in a cubic close-packed configuration, k (e.g., 
layer 'B' in the sequence ...A3G...). Some substances, 
in addition to existing in pure cubic or pure hexgigonal 
close packed structures, also exist in structures result- 
ing from a repeating sequence of layers in v7hich a layer 
within each repeat unit is in either 'h' or 'k' configu- 
ration (i.g., the sequence . . .A3CBA3C3 . . . in which a 
repeat unit consists of k layers in the configuration hldak) 
The different crystal stmcture modifications thus 
obtained are known as polytypes. Unit cells of polyt 3 np 6 s 
of a substance have the same a and b dimensions but differ 
in the c^ dimension. The c^ dimension of the unit cells 
depends upon the number of layers in sach repeat unit and 
examples where the c^ dimension of a unit cell is as Isirge 
as a fev; thousand Angstroms are known. Typical substances 
exhibiting polytypism airs SiC, Cdig > ZnS and TaS^ • Cadmium 
Iodide alone has as many as two hundred and ten known 
polytypic modifications and. a polytype of SiC with a 
repeat sequence of ^680 layers (£ dimension 12000 .X.) 

is also reported. Polytypism in various substances has 

1-3 

been reviewed by Verma and Others. 
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Polytypisa can bs considered as a one-dimsnsionaJt 
phenomenon since the only difference between any two poly- 
typic modifications of a substance lies in the sequence 
in which the layers are stacked along the c_-axis . This 
appears difficult to conprehend since truly one -dimensional 
systems noraally do not have ordered phases at non-zero 
tenperaturs . 

A polytype is completely characterized by its ABC 
layer sequence. For many polytypes, the full ABC layer 
sequence happens to be rather long and hence it is often 
inconvenient to characterize a polyt 3 pe by giving the 
full A3G sequence . Pauling, ¥yckofx and Jagodzinski 
evolved a not at ion by using whida the length of the layer 
sequence that need be specified is reduced to l/2 or l/3 
(compared to the ABC sequence) depending upon whether the 
poly type belongs to a hexagonal system or a rhombohedral 
system, respectively. 3n this system, each layer is 
characterized by 'h* or »k* depending upon whether the 
layer is in a hexagonal close -packed configuration or a 
cubic close-packed configuration respectively. The hk 
sequence now will characterize the polyt 3 pe. Another 
common notation is due to Zhdanov.^ In this notation a 
number is used to characterize a polytype. The digits 
of the number from the left, give the number of sucoes- 
sive cyclic and anti-cyciic changes (A*>B-»<3 and A:+G-^, 
respectively) in the layer orientation. 
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Polytypic aodifications of a substance form under 
almost identical conditions of growth. Often, a substance 
shows different polytypic modifications in different parts 
of a single specimen. Frequently, polytypic modifications 
show grovjth spirals on the faces of their crystals. Many 
polytypic substances exhibit random disoivier of layers 
to varying extents. Interpolytypic phase transitions 
have been obsejrved in some cases, though rarely. 

The periodicity of layer stackings in polytypes 
are often muoh larger than the range of commonly known 
atomic forces and this renders the explanation of this 
phenomenon difficult. Most polytypes are known to be 
quite stable and temperature or pressure does not 
apparently atffect their relative stabilities. However, 
teE^jerature in some cases has the affect of increasing 
the randomness in layer stackings. The number of poly- 
tsrpic phases existing fOr a substance, often far exceeds 
the maximum number of phases that can exist for a single 
con^onent system according to Gibbs' phase rule. . 

¥e mentioned eao^lier that polytypism can be consi- 
dered as a one -dimensional phenomenon. Jxi fact, we can 
consider polytypes as different ordered states of a spin- 
half Ising cXiain. The two states of the spin can be 
taken to represent the two lowest energy configurations 
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of a lajsr in a polytyps viz., the cubic aisd hexagonal 
configurations. Ail other configurations of layers 
(like AAB, A3B etc.) are of very high energies and are 
not observed exp eriinent ally , even at high tesiperatures . 
Hence there should be one to one correspondence between 
the thernodynastic properties of a polytype and the 
thermodynanic properties of an appropriate Xsing chain. 

To have ordered arrangement of spins in an Ising 
chain, it 'is necessary tc have an infinite -range inter- 
action aziong spins (see section 1, Chapter I). The fona 
of infinit s-rang3 interaction most suitable for Konte 
Carlo simulations is the equivalent -neighbour interaction. 

The contribution to internal energy from this form of 

2 

infinite-range interaction is given by ^ • This 

form of infinite -range interaction is to some extent 
justified, since such forces do exist in Nature . 

Priedcl ' has shown that the major component of elastic 
interaction between atoms of different sizes, in an 
elastic medium, has the following properties: (i) the 
interaction between li ke defects is necessarily attrac- 
tive and (ii) this force is the only interaction known, 
which is truly infinitely long-range. The contribution 
to the elastic energy is given by -cCx^-x^) .where x^ 
and are mole fractions of the atoms of two different 
sizes and C is a constant. We see that this form ^ 



Identical to the contribution coming from an equivalent - 
neighbour interaction tarn: in the magnetic Ising chain* 

IS tte postulate that the layers in the hexagonal and 
cubic coni' igurations have different thiolmessos* , we 
would hava justified the use of the equivalent -ne labour 
interaction term in the Ising chain that describes 
polytsrpism. 

If ws should have different spin orderings in 
the Ising chain in order that the chain describes poly- 
ts^isn it would be necessary to add a competing short- 
range interaction term, to the equivalent -neighbour 
interaction term. Since the 'pheunann-H^ye Ising c^ain 
described in Chapter I exhibits different spin orderings, 
we considered it appropriate to employ this for the study 
of polytypisn as well, This^mann-iij^ye Ising chain shows 
coexistence of as many as five phases for some inter- 
action strengths . This property of the Ising chain can 
then be extended to polytypes to explain the coexistence 
of several polytsrpic phases under identical ccmdltions 
of temperature and pressure, without violation of Gibbs* 
phase rule. 


*It is known, for example, that ionic crystals of the 
formula AB crystallize in cubic close-packed structures 
and covalent crystals crystallize in hexagonal close - 
packed structures, ° Sines ionic and covalent radii of 
an at om are different , our postulate is not entirely 
without meaning . 
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Wiiile ws could justify the equivalent -neighbour 
interaction term to some extent it is not easy to visualize 
the origin of short-range interaction texmis in polytypes. 
The strengths of* different short-range interactions ( J and 
K of 11*3^1) in polytypes could possibly be related to 
the growth conditions. This would imply that corres- 
ponding tc a given set of growth condition s, a particular 
set of polytypes would be formed. 

In the present study, our primary interest was in 
the Monte Carlo simulation of polyt 5 rpes by analogy to our 
studies on the Ising chain reported in Chapter I. Th# 
mechanlsh for going from one system of the ensemble to 
another, in tlis case of the Xsing chain, was a simple 
spin-flip mechanism. This mechanism loses meaning in 
the context of polytypes. To arrive at therm-odynamic 
states of a polytype, from given :toitial state, wo have 
■feo devise different mechanisms for layer reetrrangameats . 
These are discussed in section 11.3, together with the 
programs. In section 11,4 we present results and discus- 
sion of our simulation studies* ¥s shall briefly survey 
theories of polytsrpism in the next section, since this 
may be relevant to understanding the present study in 
the context of the existing Icnowledge in this field. 



IX .2 « A Srisf Revis-w ot tlis pries of FolytypisQ 

Sarly theories of polytypist: atteiapted to explain 
the phenotaenon by correlating aspects such as ispurity 
content, rats of crystallization and orientations of 
neighbouring crystals during growth, with the observed 
polytypic nodifioations . Zt was seen that the range 

of irripurity content in 6 e, 15^ SiC polytypes 

differed distinctly. Similarly, evidence fcr correlation 
between the polytyps aodification forned and the rate of 
crystallization, was found for sons polytypic sodifica- 
tions of Cdl^ . Attanpts to attribute a definite ramgs 
of growth t ssperaturs for each polytyps also met with 
sons success. Howsvsr, as nore and laore polytypic 
modifications of substances were discovered, there 
feorrslat ions 1330223 increasingly difficult , prompt ing 
& search for a more basic theory of polytypisia, 

8 

Ransdell and Sohn suggested that different poly- 
typic modifications were formed by stacking basic poly- 
meric unit (s) of the substanbe. They proposed seven 
basic polymeric units for SiC, with Zhdanov’ symbols 
33 » 32, 23 , 22, 3 k, 43 and 44. These basic polymeric 
units were arrived at from a knowledge of the stacking 
sequenoss in the known polytypic modifications of SiG. 

It was assumed that each polymer had a range of temperature 
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over -whicii it waa stable. These stability ranges soiss- 
tisjes overlapped, to give rise to two difTerent polynerio 
tinits in different eoneentrations, resulting in the 
foraetion of long -ordered polytypes. This theory was 
useful in solving the crystal structures of newer poly- 
types, However, there is no experimental evidence for 
the e^iclstencs of such polyTaeric units. Also it does not 
attempt to explain the origin of the oi*d©red stacking of 
different polymeric vaixts in long -ordered polytypes. 

9 

Verma, while studying SxC crystals, to render 
support to Frank’s screw -dislocation theory of crystal 
growth, observed growth spirals on the crystal faces. 

The measured step height of the growth spiral of 6 h SiC 
was found to be equal to the unit cell c^-dinension. 

Tnxs led Frank to suggest that polytypism in SiC is 
brought about by the spiral growth of crystals around 
screw dislocations of different Burgers vectors. Frank 
proposed that SiC crystals initially grow into thin plates. 
These plates then becosie self-stressed through nonmiifona 
distribution of impurities or thermal stresses and buckle, 
resulting in the fonnation of a screw dislocation, exposing 
a ladge on the surface. Subsequent growth of the crystal 
will have a structure corresponding to tfcat of the ledge 
and will repeat with a period equal to the pitch of the 
Screw dislocations, with Burgers vectors ifhbse 


s crew . 
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na^itudes are not integral raultlplss of the height of 

the parent unit cell, ■will thus give rise to polytypic 

modifications different from the parent crystal. Support 

for Frank's screw -dislocation theory of polytypism has 

come from the obsearvation and measurements on scre-w 

dislocations in polytypes. Direct correlation between 

step -heights of gro-wth spirals and unit cell heights of 

the corresponding polytypes have bean demonstrated in 

many cases. Another important experimental support of 

this theory comes from the observation of . syntactic 

CO ala seance (two different polytypes growing on the same 

specimen) . S3mtactic coalescence is expected since this 

theory provides for the growth of one structure from a 

screw dislocation in another. Different mechanisms have 

been proposed for the foarmation of screw dislocations 

(with Burgero vectors who se magnitudes are non-integral 

multiples of the parent unit cell hsi-^t) in different 
1112 

substances, * e.g,, impurity mechanism is favoxrred 

for the formation of such a screw dislocation in ZnS. 

13 l4 

Jagodzinski ’ raised the following objection^ 
to Frank* s screw -dislocation theory; (i) In close-packed 
structtires, edge dislocations are enargatically favoured 
over screw dislocations and the former would coup let ely 
destroy any order created by the latter. (ii) The high 
strain energy required for creating a screw dislocation 
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can .comm only wlian tbe crystal has grown to a eonsidarable 
volume by which time it would have already settled down 
to a oartain structure* (iii) While the screw disloca- 
tion s^chaaisiii requires the resulting structure to be 
perfectly ordered exar:5>les of atructures with super- 
posed random disordapr are not tmcommon* (iv) The screw 
dislocation mechanism dees not account for the theimial 
stability of polytypic etruotures. (v) A large number 
of polytypss are known, where correlation between step- 
height of growth spiral ajid the unit cell height, is 
not found* Objections (ii) and (v) have been refuted by 
Vand and Hanoka on the following basis* The greatly 
strained regions in a real crystal are not occupied by 
atoms, thus giving rise to a cylindrical hole (as 
observed in Gdl^) sround the disj^ocation* Absence of 
syntactic coalescence was not oonfixmed in the crystals 
(under (v)) the obsemred lack of correlation may 

'W&JJL be due to this* 

Jagodzinski proposed a thermodynamic theory of 
polytypism* He assumed that the internal energy contri- 
bution to free energy differences, among various polytypes , 
is mgligible * Thus the stability of different polytypes 
should be determined from the entropy contributions. 

While configurational entropy favours a disordered 
arrangement of layers, ^he vibrational entropy of the 
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system favours ordering. The total entropy, estimated by 
Jagodzinski, as a function of the degree of disorder, o , 
showed two maxima, one at o =s 0 (corresponding to a 
stable state) and the other at a -x. 0.12 (corresponding 
to a metastable state). This ic^j lies formation of 
perfectly ordea^ed polytypes as well as polytypes with 
a certain degree of randomness. Sixtytwo of the one 
hvindred and fifty crystals of SiC (the hundred and fifty 
crystals were picked at random) showed existence of 
one^dimensional disorder, rendering support to the thermo- 
dynamic theory. However, this theory is not without 
drawbacks. Long -ordered polytypes without faults have 
been observed, while it is known that the contribution 
to vibrational entropy decreases with increasing perio- 
dicity. The existence of maxima on the entropy versus 
disorder curve is based on purely qualitative arguments* 

*1 tf 

Sclineer . proposed a H^ory of poly typism> similar 
to 3 ragg-¥illiams order-disorder theory of alloys* The 
two possible layer configurations and are repre- 

sented by the two states of a half-integral Ising spin* 
Assuming that theire is an energy difference between the 
layer in *h^ and *li:* configurations , Schneer obtained an 
equilibrium ratio (d) of layers in ‘hr configuration to 
layers in configuration* Once this ratio is deter- 

mined, y interaction between the layers in different confi-* 
gurations is switched on and the arrangement of layers for 
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wh-icls the dLnteraction energy is a ainiasua is the equili- 
brium arrangement of the layers at the given teaporattire. 
Sohneer proposed that different polyt3rpes were obtained 
from a quenching of such equilibrium arrangements of 
layers. Thus polytypes, according to this theory, are 
metastable states formed on quenching the layer system, 
when the layer system is undergoing a second order phase 
transition from a cubic to a hexagonal ( 4 h) state. However, 
there is little experitaental evidence to support this 
theory. It has not been possible to associate definite 
polytypes with specific temperatures and inteipolsrtyplc 
transformations are not very common. Tbs high temperature 
modificaticai of SiC and ZnS arc not characterized by 
D = 1/2, contrary to the prediction of the theory. 

12 

Psibst * proposed a theory of polytypism according 
to which different polytypes are formed under different 
conditions of crystal growth, Ke postulated that thermal 
vibrations, under specified growth conditions, can, bring 
about a periodic introduction of stacking faults dLn the 
basic structure, resulting in the generation of polytypes. 
The observation of spiral growth, according to him, is 
associated with circularly polarized transverse elastic 
waves which may exist at the time of growth. BovovoTf 
his theory is very qualitative in nattire and offers little 
scope for experimental verification. 
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¥o see that none of the above theories (barring 
Frank’s screv -dislocation theory to some extent) explains 
all the basic features of polytypism viz,, (i) existence 
of different types of tmldimensional order while one- 
dimensional systems, usually, should not show any long- 
range ordering , (ii) seeming violation of Gribbs ' phase 
x*ule due to the coexistence of various thermodynamically 
stable phases under identical conditions, (iii) ssmtactic 
ooclsscence and (iv) existence of differing amounts of 
disorder in many polytypes. Sine© Theumann-H^ye Ising 
chain shows the above behaviour, we have attempted to 
simulate polyt 3 rpes based on the interactions of this 
Ising chain, ¥e shall discuss the ^edSkaaisws for layer 
rearrangements and also the simulation procedure in 
the next section. 

II. 3 j^ograms md Mech^anisms for lajrer Rearra^e^nts 

Monte Carlo simulation is carried out on a chain 
of one thousand layers. The calculational procedure 
involves generation of a random initial chain. Changes 
in configurations of layers are then carried out to bring 
the chain to equilibrixim at the given temperattire . Once 
the chain is brought to equilibrium, estlsiates of thermo- 
dynamic properties are carried out . At the end of the 
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calculations an output of th.e Tnlcrostates is obtained. . 

This is further processed to check for the type of 
ordering existing din the system at the giveoi temperature. 

The input parameters are, values of short-range and long- 
range interaction strengths, teaqperature and a parameter, 

A} A gives the energy separatim (in the absence of any 
interaction) between layers in 'h' or ’k* configuration 
and all other high energy configurations (corresponding 
to AiiB, AiB3 etc.) whose energies are all taken to bo 
zero. The value of A should be large enough to avoid 
o^urx’ence of higher energy states at equilibrium even 
at high tei::^6ratures. The initial state of the chain 
is always taken to be raiKlom assuming that the polyt 3 rpas 
are grown from melt.* 

The disoidered chain is created using random numbers 
between O and 1. Numbers '1», *2* and ’3’ used to 

designate A, B and C layers respectively. ¥e begin with 
the first site and call a random number between 0 and 1. 

If the random number belongs to the interval 0-1/3. *1' 
is assigned to the site. If the randoni number belongs to 

* One can, in principle, study the formation of polytypes 
under different growth conditions by either starting with 
a different initial state or by building a chain ^ 

at equilibrium such that the chain approaches equilibrium 
at each instant. The former would, correspond to growth 
of polytypes from an already existing polytypic phase , 
and the latter to a simulation of growth of polytypes 
from solution or vapour deposition. 
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the interval 1/3-2/3» *2* is assigned to the site, 
othexvise ’3* is assigned to the site. This process is 
continued till a nuaber is assigned to each site. The 
initial hk sequence is extracted froa the random A.3C 
sequence. The configuration 'h' is designated by 'I* 
and the configuration ‘k» by •2*. All higher energy 
configurations are designated by 'C . 


The initial energy is calculated using the hk 
sequence. The configuration *h' corresponds to a down- 
spin (s^ = -1) state and the configuration 'k* corres- 
ponds to an up-spin = +l) state while all other 

configtjrations correspond to = 0. The initial energy 
is calculated from the following expression; 


1000 1000 

I S^S 


1000 „ 1 
+ K I S.S. _ - AlT ^ A I 

i_1 ^ ^+2 i*l 


(11.3.1) 


Here, H is the separation between ’h* and ‘k* configura- 
tions (equivalent of the magnetic field) in the absence 
of interactions, & is the parameter already defined, 

J and K are short-range interaction paoramsters, A is the 
long-range interaction parameter and M the difference 
between the number of layers in ’h‘ state and the number 
of layers in ’k* state. A cyclic boundary condition is^^^^^ ^ 

used in our calculations . ^ ' 
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To bring the chain to equilibrium at a given tesgj- 
erature, we have to rearrange the layers. There is no 
straight-forward mechanism, like the spin-flip mechanise 
in the magnetic Ising chains, for bringing about the 
reairrangsment of layers. Although we are tempted to 
deal with the hk sequence and use the mechanism in which 
states h and k flip to each other (in analogy to the 
spin-flip in magnetic Ising chains) , a closer scrutiny 
points to the fact that such flips eire not physically 
meaningful. This is because a flip in the configuration 
of a layer changes the configuration of the neighbouring 
layers in ouch a way as to result in high energy confi- 
gurations (e.g., flipping the state of the underlined 
layei^ in tbs sequence A3C AC3 ABC leads to ABC ABB A3C) . 

It would, therefore, be best to deal with the ABC sequence 
directly and evolve physically reasonable mechan-isms 
for changing the configurations of layers. Our aim is 
to demonstrate that there exists a physically meaningful 
mechanism for rearranging the layers in a poly+ype , which 
is as efficient as the spin-flip mechanism of magnetic 
Ising chains. Once this is done, we can transfer all our 
results of the magnetic Ising chains of Chapter I to the 
case of polytypes • 

'4 . 

*Tiis results of magnetic Ising cbain ca 3 anot be transferred 
to polytypes, if no physically Tneaningf ul neclianism, as 
efficient as the spin-flip mechariis® of magnetic Ising 
chains, exists because this v/ou Id imply that changes in 
configuration allowed in a magnetic Ising chain are no 
longer allowed in polyt 3 rpes* There is really no one to 
on 0 correspondence between the states of the poly type 
chain and the states of magnetic Ising chain, over which 
thermodynamic averaging is carried out. 
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Having known the sechanisn by which changes in 
ccniiguration of the layers can be effected, we can bring 
the chain to equilibriuc at the given tanperature .by the 
probability prescription given in sectior 4 of Chapter I. 
Thsmodynanic properties of the chain are sstinsated, as 
in the casa of nagnetic Isi3ig chains , by averaging over 
a large nunber of states of the chain at equilibrium. 

¥6 have examinsd thras different mechanisms for 
rearranging the layers in tha polytypic chain. The 
simplest of thsoa, was the single layer mechanism. In 
this mechanism, a layer can change from its initial state 
(given by A, 3 or C) to any of the two other states, 
s.g., ABOARD A3C1^ . This is physically meaningful 

and ie;^li5 s that the layer in question changes its 

‘ Vs:' ■ ' 

state by shearing with respect to its neighbours. 

However, this macha-isa failed to give different statss 
of the chain, starting from a given stats, since any 
two low energy states of the ohaiipL were invariably 
connsctsd only through high energy stp-ies. This made 
low energy states inaccessible to one another. The 
low tsE^sraturs stats of the chain for negative H and . 
zero J, K and A (eq. 11.3*1) contained ’h’ configura- 
tions far in excess of tl^ number expected at equili- 
brium, (in one case the expected number was aroiand 
while the observed number was 'wl6^.) The chain seemed 
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to bs loclied in a particular state, -fax’ fron equilibriun. 
Hanes this nschanissi was discarded. 

The second nechanisn that we experisientsd with 
was tha layer transposition •mechanisai due to Jagodzinski. 

In this isGchanisBi, two adjacent layers in the chain inter- 
change their ABC configuration, e.g., A3G AC 3 ABC 
A3C ABC ABC . This mechanisn is also physically reasonable 
since it involves a shearing of the two layers in question. 
This nechanisn was proposed by Jagodzinski to explain 
formation of different polytypic nod ifica.t ions froa an 
ordered phase. ~J°. carried out a siaulation using this 
nechanism for J >0 and H, K and A equal to zero (Sq.II.3.l) 
¥e found that the nunber of like neighbours (hh or kk) 
was far larger than the nunber expected at equilibriui:. 

A closer look at the aiorostates showed that we were 
led to configurations from which we cannot come out, 
except via high energy states (e.g., for J > O & A, K & E=C, 
ws observed states like ...ACBCA^SO, which can never go 
to ACBCAC3CA by layer transposition mechanism, though the 
latter is a lower energy state). Jagodzinski* s mechanism 
thus failed to yield a polytypic modification, starting 
from a disordered state. A combination of the single 
layer mechanism and the layer transposition mechanism also 
similarly failed. To overcome these difficulties. 


we 
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brought the higbsr energy statss slightly lower emd 
carried out a gradual cooling experi^sent* . ¥e were 
once again raced with the sacse prcolen viz., locking of* 
the chain in non-equilibrium configurations unless the 
higher energy states were brought sufficiently low so as 
to be accessible at low temperatures, a situation in 
contradiction to experinaital observations. Using these 
Eschanisns, in the case of A, K and J equal to zero and 
li less than zero, the observed number of 'h' states v/as 
3«5'7& while the expected number was 'v-O^. 

A successful mechanism of layer rserrangement 
should be capable of connecting two low energy configura- 
tions directly (i.e. , without an intermediate high energy 
configuration) , This can be achieved by a general two- 
layer mechanism. In this mechanism, two adjacent layers 
Eire allowed to cbemge their layer configurations arbitra- 
rily {e.g., A3CAGB -» A3C3C3 or ABABG3 and other high 
energy configurations) . Changes to high energy configu- 
ration are highly iaproba.bls and changes from one low 
energy configuration to another low energy configuration 
alone occur at the temperatures of interest . This ^ 
mechanism is physically meaningful since such changes can 

Eere , we start from an arbitrary high temperature state 
and lower the temperature in convenient steps, till the 
desired low temperature is reached. At each te 5 %>sratur©, 
we take the chain through a sufficiently large number 
of Monte Carlo cycles so that the chain comes to 
equilibrium at each of these temperatures. 
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be brought about by shearing the layers such that the 
atcr.s in the layer move , always keeping to the hollow 
space formed by the atoms in the neighbouring layers 
(during the whole process these atoms in adjacent 
layers never come one above the other while going from 
one low energy configuration to another) . This mechanism 
gave nearly equilibrium number of hh and kk configura- 
tions, for J>0 aoid K, £ and A all zero (eq. 
in a reasonable number of Monte Carlo cycles. Using 
this mechanism, the computer time taken (on IH'I 7044/l4ol 
system) for carrying out 120,000 Monte Carlo cycles was 
about 2C minutes. A display of a section of the chain 
at different stages in its approach to equilibrium via 
the double layer mechanism is given on the next page 
(Fig. II. 1) . 

It is difficult to see the existing order in the 
chain, in any given stats, just by looking at the chain. 

The order parameters in the form of correlation functions 

■ 'ife 

are obtained by an averaging of the same over a large 
number of states of the chain and hence do not give the 
ordering of just a single state of the chain. To search 
for order in a given state of the chain, we have employed 
a display program. In the display program, we start 
with a lanit of a given number of layers, at the beginning 

This inplies that a high activation energy barrier' that 
was separating any two low energy states in our calcula- 
tions has been overcome . Thus, we can achieve equilibrium 
configurations in much fewer Konte Carlo cycles. Ihysically, 
this means that equilibrium is achieved in shorter times 
and the chamges in configurations qre nc longer kinetically 
controlled but are controlled by thermodynamic probabilities. 



Figure II. 1. Growth of 4H Polytype Using the General Ivo Layer Mechanism: (a) After 4'';000 Fonte Carlo 
Cycles (b) After 95000 Monte Carlo Cycles (c) After 145000 Monte Carlo Cycles (d) ifter 195000 Monte 
Carlo Cycles. 
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of' tb' chain and s-ie if* it ns-^eats itssl-f. I:? it does, 
then we continue till there is a iniciEatch. (stacking 
fault). At the misBiatch a blank, (or a dot) is stored 
and ws proceed with the comparison using the same number 
of layers, beginning immediately after the mismatch- If 
there is a taismatch, even before the set of layers has 
repeated once, ws store a blank (or a dot) in the first 
layer of th.s initial set and start afresh with a set of 
layers (containing the given number of layers) immediately 
after the mismatch and proceed in the sane fashion till 
we rsa^ch the end of the chain. ¥3 can search for order 
of a.:ay periodicity in a polytyps chain (or a magnetic 
Ising chain) by employing this program. In our displays, 
ws have searched for order in hk seguenoe and have later 
c onvsr t sd them into the ABC seausnes . 

II. 4 Results and Discussior . 

The Monts Carlo simulation of poly types reported 
here is based on the Theumarn-H/iiys Ising chain. The 
interaction strengths are so chosen as to cover all 
regions of Theumann-H^^ye phase diagram ( Fi g . l.l), with 
one set of interaction strengths per region (actual 
values are the same as in Table 1.1 ). ¥© find that the 

actual numerical value of A (defined in the previo^^ 
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sectioR) does not alter the results of* our calculations, 
as long as its absolute value is muci" larger than the 
absoluts val’ss at interaction strengths. In each region 
o:r the Theusann-Ii/ye phase diagram, we have carried out 
a Monte Carlo simulation for a single teiaperature , this 
temperature corresponding to the lowest tanpera.ture in 
our calculations on the magnetic Zsing chain (Chapter l) . 
In ail our calculations, H (Sqn. 11.3-1) is taken to be 
sero; mn-zero li would imply that, in the absence of 
interactions, the layers in the hexagonal and cubic 
configuration s differ in stability. 

In region I, the polytypic modification obtained 
was 3C (cubic close-packed structure). The 2H nodifica- 
tron (nexagonal oloss-packsd structure) is also equally 
likely sines for zero H, there is no preference between 
'h* and 'k* configurations. 3y starting with a different 
random chain, we could well have obtained a 2H polytype. 
¥e would have expected these polytypic modifications, 
from our results on the magnetic Ising chain, for this 
region. A display of the actual polytype is given in 
Fig. II. 2. Ye should also expect this polytype to shov; 
a first order phase transition. The high temperature 

* 2H and 3G are tbs two extreme polytypic sodificatipna 
of a substance. m;. 



Pigure II. 2. I/isplay of lolytype Formed in Hegion 
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piiasG would ha-wB a sliort-range ordering o£ the type hhlch 
wiiils the long-rangs order still remains the same (Fig. 1.7) 

In region II, the low tsmperaturv; ordering ox the 
Isyers corresponds to a 12E polytyps (ABAC3C3ACAC3) . This 
is whet ws should expect irom our studies of the magnetic 
Ising chain, since spin configuration corresponds to 

hh kk or kh hh layer configuration.! both these lead to 
tii 3 sans 12R polytype . A display of this polytyps, 
obtained from our Honte Carlo simulation, is given in 
Fig. II. 3 . In this region, the magnetic Ising chain 
showed two phase transitions and the same wotsld be expected 
of the pclyt 3 rpe chain, "ife showed that , for the magnetic 
Ising chain, the short-range oixier remains essentis.liy 
the same , while the long-range order changes, as the 
chain undergoes a phase transition. The sane should be 
true in the case of the polytyps as well. From Fig. I.13» 
we see that this transition implies a change in long- 
range ordering from 12R to 9R (.ABC3CACAB) in the case 
of the polytype . 

In region III, the polytype does not have a single 
ordering of layers throughout the chain . In some parts 
of the chain, the polytypic modificat ion is 12H while in 
other parts i the polytypic 5iodificati<m is either 
(3AC3CA.) or 9S (ABC^^ACAB) . Presence of different 
polyt5rpic modifications in the same specimen has indeed 



Figure II. 3- Displaj of j?olytype Formed in Ee.: ion II 
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'bsen obsei^sd exp eriuen tally ixi rr’any systes:c . 
results on nagnstic Ising chain, we should expect to sas 
several second -order phase transitions in this region 
(though W6 are not in a position to coisnient on the 
changes in ordering) . 

The polytype obtained jfron region X7 belongs to 
the 9S nociricat ion . Onca agein, like 3.11 region I, we 
could have instead obtained a 6 h polytype • The 9 ^ 

6li polytypes correspond to the ++4- and spin configura- 

tions, respectively. Display of* the polytype foriaed in 
this region is given in Fig. II. 1. The polytype f*omsd 
in region T could be considered as a 9 H modification 
with a Isjrge number of stacking faults. Similarly, the 
polytyps from region YI was a kU (A3A.G) nodificat ion 
77ith a large number of stacking faults. Experimentally, 
such polytypes with a large amount of stacking faults 
are knoi-^n to exist. Polytypes in these regions (IV, V 
and VT) arc again expected to undergo several phase 
transitions. An. increase in temperature will have the 
effect of increasing the concentration of stacking faults 
(as can be seen from the correlation functions in Chapter l) 
in conformity with experimental observations on polytypes. 

Tfe extended the range of short-range interactions 
Lip to fifth neighbour, in an attempt to generate longer 
psriod polytypes in cur simulations. In this atteo^jt . 



Figure 11.4. Display of Folytype in Region IV, in Field 
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7je experimented with different Toms of short -range inter 
action potentiels. These are shown in Fig. 11,5. Sitmie 
tions with interactions of the form ’a' and *b' (given 
in the figure) resulted in 4H polytSTPSs with varying 
amounts of disorder. It, therefore, appears that short- 
range interaction siiailar to these forms (potentials 
gradually decreasing with distance) are not responsible 
for formation of long -ordered polytypes. Simulation 
-rfith the short-raiige interaction potential given by 
Fig. II. 5c resulted in the formation of 12E polytypic 
modification (hhiikldk) . A display of this polytype is 
given in Fig. 11.6. Contrary to our expectations, an 
increase ir. the range of short-range interactions did 
not result in the formation of polytypes with longer 
periodicities (compared to tbs Theumann-H/jiye chain). 

The longest polytypes found by us wei*e the 12R and 12H 
polytypes corresponding tc a ^-dimension of about 30 
This in itself is quite encouraging considering that we 
have not rea.lly focussed on carrying out calculations 
with various possible combinations of interaction para- 
meters. It may be worthwhile to examine whether there 
sire other forms of short-range interactions (together 
with equivalent-neighbour intera.ction) which give rise 
to poiytypes with very large periodicities cocpared to 
those encountered in the present study. 




Figure II. 6, Polytype Fomed From the Short-range Interaction of Fig. Ifel-Sc 
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11.5 Suini:ia^ 

Since polytsrpxc mocixfications of a substance differ 
from one another in the stacking sequence of layers in one 
direction [ c^) , we can treat then essentially as one- 
dinensional syctens . A polytype can, in fact, be consi- 
dered as a spin -half Ising chain since only two possible 
oonfigurat ions for a layer exist, viz., the hexagonal 
and the cubic. Existence of polytypes seems to contradxct 
Gibbc^ phase rule because polytypes far in excess of the 
naximum number of phases allowed by the phase rule are 
Icnown to be stable under identical conditions of tempera- 
ture and pressure . Fuxiihezmsore , in one -dimensional systems, 
unless there is an inf ini te -range interaction, there does 
not exist any long-range order above absolute zero. Thus, 
polyt3rpes seem to have a behaviour similar to tnat 
exhibited by the Thcuranr-H/eiyB Ising chain studied in 
Chapter Z. The long-range interaction term used in th© 
Thsumem-H^eiy© Ising chain is justified for polytypes , if 
we assume that the layers in the two configurations *h’ 
and have different thiclnisssss. Although we do not 

have any real justification for using the short -rajxg© 
interaction terms of Chapter I , we feel tnat this may be 
related to growth conditions . 
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To carry out a Mo.”ts Carlo sinulatior of polytypoo , 

*7© nact have a neclianism for rsarrarging the layers ±n 
analogy to the spin-flip nochanissi used in aagnetic Ising 
chains. To extsnd our results of aiagnetic Icing chains 
to polytypss, we should arrive at a physically meaningful 
nschainisH capable of comacting any to/o states of a poly- 
type, with a probability which is about the same as the 
probability with which spin-flip nechanisia connects the 
two corresponding states of a magnstic Ising chain. 
have studisd three different layer rearrangement nechanismo, 
viz . , (i) singl;. layer mechanisms (e .g. , -*-A3G&BC) > 

(ii) layer transposition mechamism (e .g . , A3CAC3 AC3A.C3) 
and (iii) general two layer mechanism {e.g., AECAC3 •» A3C3CB 
or ABABCB) . ':Ie find that the general two layer mechanism 
satisfies the criteria set forth viz., it is physically 
meaningful and ergodic. 

’fe have carried out a direct simulation of polytypes 
for one value of temperature and for interaction strengths 
corresponding to different regions of Theuriann-Ejitye phase 
diagram (5'ig . 1.1 ) . ¥e obssirve different pol 3 rt 3 rp 6 S in 
different regions as expected. The various polytypes 

■ ".A 

observed are kR, 6l-I, 9R end 12R, in addition to the basic 
2H and 3C polytypss. In region III, the polytypic chain 
shovjsd 12?. ordsring in some sections (of the chain) and 9^^ 
or 6 h in some others, an observation supported by experi- 
ments (syntactic coalsscencs ) . To effectively generate 
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different poly types on diffsrer.t parte of a chain, 
should deal with longer chains and take then through 
larger number of Konts Carlo cycles. These are expected 
to undergo phase transitions by analogy with the nagnetio 
Ising chain. Associated with those phase transitions, 
we should observe changes in long-range oi^er. !Je also 
siriulated poly types with different forms of short-range 
interaction, extended upto fifth neighbour, in an attempt 
to genarsts polytypes with longer periodicities. This 
was not very successful arid the longest polytyps observed 
in this study was a 12H polytype . 

The origin of short-range interaction terms used 
in the present study is still to bs understood. However, 
an understanding of the mechanism of layer rearrangements 
together with Monte Carlo simulations provides a nsans of 
directly verifying the formation of polytypjss from a 
proposed interaction model, 'fe can also use this method 
to simulate different sets of growth conditions and study 
the formation of polytypes under any given set of growth 
conditions. In the present stage of our study, we have 
not concerned ourselves with the formation of growth 
Spirals on the crystal faces of most of the polytypes. 
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CEAPISR III 


models FOR LO¥ SPIN -HIGH SPST TRAlISITICae 22T SOUmS 


111,1 Introduction 

IJlien a transition metal ion is placed in a crystal 
field, tile sliift in energies of the various terms, with 

changes in crystal field strength is givm by the Tanabe- 

1-3 f h 5 6 

Sugano diagi*ams. For some configurations (d , d , d 

7 

and d'), the diagram predicts a crossover in the ground 
state from a spin state with a lower multiplicity to a spin 
state with a higher multiplicity, with change in crystal 
field strength. In systems where the crystal field strength 
lies close to the crossover points, a t he rnKil equilibrium 
may be established between the low-spin states and the high- 
spin states. Such a spin-state equilibrium in a system gives 
rise to inteiresting magnetic behavioar. Existence of a spin- 
state equilibrium was first established in the iron (iTl) 
oo!35)lsx , Fe (dimethyl dithiocarbamate)^ .^ 
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Recently, sosie dooardination coispotmds of iron (IX) 

f 6 

(d configuration) hare been studied in detail in the solid 

4-7 

state. These cooplexes seem to tmdergo gradual or abrupt 

changes in the spin-state populations, often accorrpanisd by 
a phase transition. The nature of this phase transition is, 
however, not well understood. Spin-state transitions have 

8 9 

also been obseived in inorganic solids such &s IhiAs and IfcP, 
COgO^,^^ LaCoO^,^^’^^ HoCoO^^^ and other rare earth ortho- 
cobaltites. Spin-state transitions in some of these ino::ganic 
solids, unlike in the case of the coordination conplexes, are 
associated with structxzral phase transitions. particularly 
interesting are the spin -state transitions in LaCoO^ and 

HoCoO^. The transitions in these two systems are associated 

1 h- 

with an unusual structural phase taransition in which the 

V i ■ ' 

high-temperature phase is of lower symmetry than the low- 
temperature phase. The low-tempeaTature phase in the case 
of these cobaltites, has a residual high-spin population. 

The spin-state transition in these cases is either gradual 
or abrupt and in the high-temperature region, the populations 
of the low- Spin and high-spin states are fotand to be nearly 
equal. The spin-state transition in the cobaltites is 
followed by an ordering transition resulting in the symmetry 
change , In the ordering region of I'faAs , I'SnP and the rare 


* 15 

PrAlo^ is another system in which the high- temperature 
phase has a lower symmetry than the low-tenpjerature phase , 
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©ajrirli oobaltites a plateau is obsearved in. tbe inverse magnetic 
susceptibility versus tempearature cuarve . 

Althougb some macroscopic approacb.es have been ©n^loyed 
to explain spin-state tjcansitions in solids with some reason- 
able success, particularly in tns case oi cooardinaticai 
compotmds, a satisiactoiy microscopic model for these phase 
■ tai^nsitions is still lacking. To explain the phase transition 

16 17 

in LaCoO^ , Bari and Savardiea?© extended Chestnut's theoary 

for first order sing let -trip let transitions in radical salts 
IS 

of TCNQ . Bari also proposed a two-sublattice model to 
eaqslain the change in symmetry following the change in spin- 
state population. However, these models are not successful 
in explaining the observed behaviour near the transitions in 
LaCoO^, HoCoO^ and related systems, A successful model for 
spin-state transition should be able to explain the following 
observations; (i) smooth as well as abmpt changes in the 
spin-stats population ratio with temperature, (ii) the occur- 
rence of thermodynamically first or second order stxT^ctural 
transitions in the system and (iii) non-zero population of 
high-spin states at low temperatures, found in some systems. 
Besides in the rane earth oobaltite systems , plateau in the 
inverse magnetic susceptibility versus temperature curve as 
well as the unusual structural transition resulting in a 
lower symmetry high-temperature phase also need be explained. 



Because of* consideoTable experinieritai interesi; in "bliis 
laboi'atory in spin-state transitions, we decided to examine 
soma models to explain such a transition. In this chapter, 
we. present results of our investigations of some new models 
for high-sp in-low -Spin transition in solids. ¥e have been 
able to show that the coupling assumed in 3ari*s model is 
incapable of producing a phase transition if wer consider the 
dynamics of the lattice. To explain residual population of 
the high-spin state at low temperatures, we have studied a 
model in which the crystal field strength is quadratically 
coupled to the lattice displacement, with the co'-pling capable 
of mixing the spin states . This model predicts a residual 
high-spin population at low tenperatures , although, like in 
the previous case , it does not exhibit any phase transition 
and only shows a gradual variation of the ratio of the two 
spin-state populations with tenperature . Since a structuaral 
transition is associated with the spin-state transition in 
systeas of our intej^est, we have obtained a free enei^y 
expemsion in teimis of displacement by employing ti» path 
integral method to evaluate the partition ftanction, (Details 
of the method are given in the Appendix.) 3n the two- sublattice 
case with linear coupling, we obseive a second order phase 
transition with a gradual change in high-spin-low-spin popula- 
tion ratio, as a function of the order pai^meter ovei* 
sublattic6 . ¥e have studied a two -sublattice model with cubic 
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coupling of the spin states to the lattice and we find that 
this systea is capable of giving a themodynamically first 
order structural transition » A two phonon model has also 
been studied by us. In this model liie spin-states are 
coupled linearly to the breathing mode and quadrat ic ally , 

•with mixing, to the optical (ion versus octahedral cage) 
mode. This model shows a phase transition and the ratio of 
Spin-State populations depends upon the. order parameter. 
]rurther it also shows a non-zero high-spin state population 
at absolute zero. 

In the next section, "we give brief review of the 
experimental situation in spin-state transition in solids . 

Ill .2 Spto-State ‘^ansitipns^ in Solids - The Experiment al 
Situati<m 

'■'W" 

Spin-state transitions in coordinaticsi coEiplexes of 
iron (ll) have been studied extensively. Sorai and Geki* 
observed heat capacity anomalies in the iron (ll) con^jiexes 
I|'e(phen)2(NCS)2l] aaid [Fe(phen)2{NCSe)23] at 176 s and 231 E 
respectively. The inverse susceptibility versus temperature 
plots show a plateau aroinid the ten^eratuires at which the 
heat capacity anomalies are obseived, .The respective effec- 
tive magnetic moments show abrupt jun^js at these temperatures. 
The infrared spectra of these compounds show intei^stihg 
variations with tes^ei^ture. The spectra of the sulphur 
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coii5>1®3c sIiow a doublet ( v_ ) centred around 20^8 csa*”^ at 

o— rt 

low teiaperatures , As tbe t enpe j^ature is increased, anotliar 
doublet centred arotmd 2110 da"^ appears. Tba intensity of 
tbe low -teaipe nature doublet decreases witb increasing 
ts 25 >e 3 ratu 37 e while the high-temperature doublet shows the 
opposite taTend, The positions of the doublets remain 
unchanged with temperature. The results are explained on 
the basis of a spin-state transition in these systeas- The 
infrared studies indicate that the energy separation between 
the two spin states remains the same at all temperatures. 

Direct evidence of spin-state transitions have been 

4 5 

obtained by Mossbauer studies. Low -spin and high-spin 

Fe(U) ions have different isomer shifts thus making possible, 
separate quantitative study of t be individual spin, states. 

The iron (in) complex [>e(4,7-(CH^)2-I’li®»)2(^^^)2^ 

fi^(Papt) 2 J (papt - 2 _( 2 -pyradilamlno)-^-( 2 -pyridyl)thia 5 solate) 

. >■ 

have been studied enploying M 6 *ssbauer spectroscopy. The 

studies reveal that [Fe{ 4 , 7 -(GH^) 2 -P^®ii) 2 (^'^^®) 2 ^ 
an abrupt low-^in-high-spin transition while [ 5 '®{Papt) 2 “| 
undexgoes a smooth spin-state transition, |i'e(papt) 2 J also 
shows a residual population of the high-spin sta.te at low 
tesperaturss . The spin-state population in []^ (4 , 7 -(CE^) 2 - 
Ph 6 n) 2 (NCS) 2 l] is known to show hysteresis behaviour with 
temperature^ (indicative of a structtiral ti'ansition) • However, 
no definitive structural transition appears to ha.ve been 
rspojrted to accompany a spin-state transition in any of the 
coordination conplexes. 
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Amca^ isimple inoiganic solids, MnAs and are tiie 

first solids in •which spin-state transitions were observed - 
MnAs eshibits "two stmctnral phase transitions -which are also 
associated -with chaises in magnetic properties. The phase 
transition around 315 K is a first order transition -with a 
change in structure from B8^ to B31 . The second phase transi- 
tion occurs around 400 K and the struc-feure changes fron B31 
to • The low-temperatiirs B8^ phase is ferrosiagnetic 

=3.^ Pg) and the high-^tenperature 38^ phase is para- 
magnetic = ^#95 Vg)» 'I‘3ie 331 phase does not exhibit 

any magnetic ordering and in this region the inverse s*uscep- 
Ifl'y ''^®rEus temperature curve has a negative slc^e. 

These results have been qualitatively explained by Goodenou^, 
Hidgley and Neuman,^ on the basis of a spin-state transition 
brought about ■when the energy pajpaneters occurring in the 
band diagrams ( of l^nAs in B8^ and B31 phases) satisfy cesrtadn 
inequalities. 

Cobalt sesquioxide, 0020^, is another system where a 

transition has been observed. The ionic i*adius 

of liigh-spin cobalt (Hl) is significantly largier than the 

ionic radius of low-spin cobalt (XXl) . 002 ©^ synthesised at 

pressure ^ thus has all the cobalt icGis in a low— spin 

state. On annealing COgO^ sit 673 K» the crystal changes 

from a hexagonal structure to a co3cundujn structure with 6 , 7 ^ 

increase in volume. Con 5 >ariso(n with X-ray patteju. of a-Fe« 0 - 

2 3 
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indicates that all the cobalt ions in the corundum structure 
are in the high-spin state. Although this transformation has 
not been reversed or fully characterized, it appears that 
this transition in COgO^ represents a first order spin-state 
transition. 

The rare earth cobaltites that have been studied so far 
can be classified into two categories. The first category in 
•which the lo'w-spin Co^”^ ions first transform to the high- 
spin state upto a temperatinrS , followed by fornsation of 2^ 
and 4'*' oxidation states of cobalt arising from charge- transfer 
be'tsi/een the t'wo spin states of trivalent cobalt and the 

■ 'O ; 

second ca.tegory where low-spin Co"^^ ions transform to the 

high-spin state upto a temperature and then exhibit a constant 

proportion of the two spin, states upto high temperatxrpes . 

11 12.24 

The cobaltites belonging to the first category are LaCoO^, * ^ 

19»24 _ „ 20,25 ^ 19,24 ^ „ 21,24 _ 

NdCoO^, EuCoO^, ’ GdCoO^ and YGoQ^. The 

13 24 

cobaltites that belong to the second category are HoQoO^* ’ 

2 1 24 22 22 2 

SrCoO^, * YbCoO^ and LuGoO^. * ^ In addition to these, 

cobe.ltites of prase odyaium and Dyspr osium have also been 
24 

studied. In Table 111,1 •we have giv«i a brief summary 
of Spin-State transition stiadies carried out bn the various^^^^ 
rare earth cobaltites (incliKiing YCoO^). Only cobaltites 
belonging to first category show thermal anc^alies (e.g., 
differential theres-l analysis curves) around the spin-state 
transitions inplying that the spin-state trans ititm is followed 
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T^le I 

Spin-State Transiticms to Kara Earthy Cobaltates 
Transition ten^ei^tures, E f^p^ 



DTA 

X-ray®^ 

Magnetic 

suscepti- 

bility}’ 

L-M ^ 

factor^ 

Center ICo^'*' ! 
shift‘d - 

f 

AH , Kcal 
mol~1 

LaCoO^^^’^^ 

680 

500 

620 

650 

650 

1.0 

0.2 or much 
greater^ 

NdCoO^^^’^^ 

600 

580 

600 

- 

4 oO 

0.8 

<0.2 

EuCoO^^®"^^ 

- 

600 

600 

600 

600 

0. 8-0.9 

<0.1 

GdCoO^^^*^^ 

530 

450 

550 , 

500 

600 

1.3 

0.1 -0.2 

YCoO^^^ 

650 

700 

520 


650 

0.8 

0.2 

„ „ ^ 13,24 
HoCoO^ ’ 

- 

325 

300 

300 

300 

1.0 <0.1 
(rises from a 
non-zero -value) 

^ 21 ,24 

Ei>CoO^ ’ 

600 

600 

600 

- 

550 

1.0 <0.2 
(rises from a 
non-zeiTo value) 

^ ^ „ 22,23 
LuCoO^ ’ 

- 

54 o 

450 

- 

- 

1.0 <:0.1 
(rises from a 
non-zero value) 

_ _ _ 

, 

.. 



▼ 

^ ^ V . ^ 

^ ^ nn*. 


From Debye -Trailer factor data* Temperature cori^spondtog to tlae 

beginning or end of plateau region to the 1/x~temperature curve j 

c d. 

from Lamb-Mossbauer factor data; Tcmperatuare wiiere tbe Mossbauer 

0 

centre shift decreases marlsedly; Relative proportion of bi^-spto 

3 + ' ' . ' ' ' ' ' 

Co at DTA transition temperature; values from DTA; ' &E seems to 

depend on method of preparation. The transition could indeed be 
t her modynaniic ally first order, particularly to LaCoO^ prepared by 
coba It i cyanide decon^ositicm. 
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by & phs.se transition. The &E values of the endotlieraic 
transitions (whan observed) are generally snail (Table III.l); 
sucii AH values could ia5>ly occurrence of second order transi- 
tions. Anotlier general feature of tbe spin -state transitions 
in tns cobaltites is tbat around tbe transition tenperature , 
tlie ratio of higb-spin population to tbe population of otber 
Spin states is unity (Table UI.l), In tbe rest of tbis 
section ws sball confine ourselves to a discussion of 
relevant studies on LaCoO^ snd HoCoO^ s-s these are representa- 
tive of tbe two typical categoiies aentionsd earlier. 


Systenatic studies on LaCoO^ were first carried out 

11 

by Sac cab and G-oocl enough and later by 3 bide , Hajoria, Rao 
12 

and Rao. These investigators found a plateau in tbe inverse 
susceptibility versus tesperatui^e curve in tbe region 
^00 K <T< 650 K. Tbe variation of inverse susceptibility 
with tecperature is fotmd to be linear in tbe regions 
4,2 K<T< 400 K and 650 K< T <1200 K, although tbe slopes in 
tbs two regions differ. Tbe unusual nagnetic behaviour 
stzggested the existence of two different spin states of 
cobalt ion in equilibrium. The DTA curve shows a broad ^^^ ^ ^ 

endotherniic pealj: around 680 K indicating a phase transition 

. ' 11 

in the system. X-ray diffraction studies shew that the 

Space group symmetry of laCoO^ changes at 65O E from R 3 G 

(below 650 E) to R 5 (above 650 K) . Below 650 S, there is 

only a single 0 -Co distance while at 650 K the 0-Co distance 



versus tsaperatura cuive brandies, indica.tang two distincrt 
0-Co distances, above 650 K. Tbe volume of' the unit cell 
does not show any significant change at this tenperattire 
though the volume versus ten 5 >erature plot shows changes in 
slope at 4oO and 650 K. Debye-Tfaller factor (from X-iay 
measurements) shows abrupt changes at 400 K and 650 E. Ih 
the ten^jerature region 4oO K < T < 650 E, short-range ordering 
of two different t 3 rpes of cobalt ions is suggested based on 
the Debye -Tfaller factor studies. 

The variation in spin-state equilibrium as a function 

12 

of tenperature was studied by 3iide et al., ecploying 

57 

Mossbauer spectroscopy taking Co diffused LaCoO^ a-s the 
Mossbauei* source and E2j^Fe(GN')^.3^0 as the absorber. It is 
assumed that the decay does not affect the spin-state of 

the ion i.e. a high-spin Co^'*’, on decay, gives a higb-spin 

3 3 3‘4' 

Fe and similarly a low-spin Co"^^ gives a low -spin Pe « 

Mossbauer spectra of LaCoO- show a doublet with isomer shifts 
0.0 +_ 0.04 naa/sec and 0,45 hh 0.o 4 aaa/sec and can be safely 
assigned to tbs two spin states of Co^^ ion. The ratio of 
the high-spin population to low-spin population has been 
obtained from intensity measurements in tbs t©::;perature legion 
4.2 E<T< 1200 E. This ratio is also calculated from the ^ 
susceptibility data and these are shown in Fig. IH. la. The 
susceptibility studies indicate that the high- spin ion popula- 
tion increases with teaiperature upto 1200 K while M2>ssbauer 
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studies show that the high-spin ion population goes through 

a jaaxijjiuia around 200 K to zero at 1200 K. This discrepancy' 

between the Ilossbauer and susceptibility results is eaplained 

on the basis oT charge transfer between the two spin states 

of cobalt (HEl) resulting in low-spin cobalt (ll) and high- 

snin cobalt (jV^ « The latter pair of ions also have an isoEter 

shift of nearly 0 sjia/sec and overlap with low-spin cobalt (Hi) 

resonance . Mossbauer studies also indicate a phase transition 

around 650 h as evidenced by (i) a aarlced decrease in the 

centre shifts of both the spin states around 650 K and (H) 

dips in Lanb -Mossbauer factor-tesiperature curves at 4oO K 

and 650 K, sinilar to those observed in the Debye-v7aller 

factor versus teaperature cuives, Sxistenoe of the low- 

spin and high-spin states of cobalt (Ul) , in equilibriina, 

27 

is also confi 2 T 2 sd fron optical studies. The energy difference 
between the low -spin and high-spin states in rare earth 
CO baltites, estimated froa various studies, range from 0*01 eV 
to 0*2 eV, 

The Spin-state equilibrium in HoGoO^ bas been sttaiied 
by Bhids, Rajoria, Heddy, ha o and Rao. The inverse suscep- 
tibility versus temperature curve of lioCoO^ shows a plateau 
in. the region 250 K<T< 300 K. This variation is linear, 
below 250 S and above 300 K, although ths slopes in these two 
regions are different. This behaviour suggssts a spin-state 
equilibrium in the system. The DTA curve of HoGoO^ does not 



show peaJc around, the plateau region and hsncs no definite 

phase transition is associated with the spin-state transi- 
tion. Detailed X-ray inArest igations have not been carried 
out' on HoCoO^ and so no definite structuial transition has 
been established, unlilce in the case of LaCoO^. However, 
Debye -Haller factor versus ten^ei^-ture curve exhibits breahs 
around 300 K and 600 K, although the X-ray unit cell volane 
varies smoothly with teaperaturs in the interval 300 K < ? < 
1200 X. 


Ilossbauer spectra of HoCoO^ show doublets, corres- 
ponding to the high-spin and low-spin states of trivalent 
cobalt, in the region 78 K < T < 800 K. From intensity 
measurements, it was found that the ratio of high-spin 
population to low -spin population is constant at '*'0 -7 » belcaw 300E 
Around 300 K, this ratio shows an abrupt jump to a value 'ul ,o 
and above 300 K, the ratio remains constant at this value. 


The Spin-state population ratio obtained from susceptibility 
measurements agrees with the Mossbauer data (Fig. III. 1b) . 
Thus, ■unlihe in the case of laCoO^, we need not postulate 
charge transfer be'tween diffei*snt spin states of cobalt (jIT) , 
in EoCoO^. Center shifts of the two spin states of HoCoO^ 
show a marlcsd decrease aroxmd 250 X and Lamb-Mbssbauer factor 
variation with tesperaturs closely follows the temperature 
variation of X-ray Debye-Haller factor. These results point 
to an ordering transit ion in HoCoO^ , similar to that observed 
in LaCoO^. 
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Temperature dependent resistivity and Sesbsck coeffi- 
cient measurement B have been carried out <ai rare earth 
cobeiltites. The variaticcis of these properties with ten^jera- 
ture reflect the changes in the nature of cobalt states found 
by Mossbauer and susceptibility measurements. The rare 
earth cobaltites are p-tj^e semiconductors at room tetqjera- 
ture and undergo a semiconductor i. metal transition at high 
temperatures ('vl200 K) , an aspect of considerable interest 
but not of direct relevance to the present discussion. In 


the next section we give a brief introduction to the models 
of Chestnut Bari and Sivardiere . 


III. 3 3arlier Models for Low-Spin-1-y.gh-Spin Transition In 
Solids 

In order to explain the first-oider transition between 

17 

singlet and triplet states of radical salts of TCa-IQ, Chestnut 
consider^ a model system consisting of N molecules, n of 
which are excited to a state of degeneracy v. The inteimal 
energy of the system was assumed to be of the form, 

I = P®o (III.3.1) 

where and are parameters* of the model and p = n/R is 
the particle density. The entropy of this system is given by, 

o/N = plnv-plnp- (1 _p) In (t -p) {ITT. 3.2) 

e is the single particle energy required to create an excita- 
tion. The term in can arise either due to particle- 
particis interaction in the excited state or due to a 
deformation of the lattice caused by the excitations . 
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Thus , tlis free energy of tlis systen can be -written as 
F 1 2 *! 

jjr = EoP + 2 £i P - g; {pin V - plnp - (l-p) In (l-pQ • 

(111.3.3) 

Free energy is sinisiized with respect to p to obtain the 
equilibrium value of P as 


p = v/ {exp [3g( p)3 + v}. 


together with the condition 


1 

K dp2 


€ 

1 


1 

^ gp (1-p) 


(111.3.^) 


( 111 . 3 . 5 ) 


¥hen the equations (ill. 3,4) and (111.3.5) are solved, for 

various values of e . e. and v >1, below a certain value of 

o 1 

(-e^/e^), p shov/s a gradual increase with tenperature while 
above this value of (-e^/e^), two equally stable minima are 
observed at a transition tenperatui^ . Above the transition 
temperature a higher value of P is praf erx*ed . The limiting 
value of p , as T is v/(v + l) . The jucp in particle 

density, from p:j to pg , a-t the transition temperature is 
characterized by a symmetry around P = O .5 . 

Bari and Givardiere^^ extended Chestnut ’ s model to 

explain low-spin-high-spin transition in solids and chose 
* 

the Hamiltonian , 

Tills Haiail'tonia.n gives exactly tiie sai^e foi* in'fesxtn.aJL 

energy^ as assiiiijed by Glaestnut (eq. IIIe3#l)# Titis can 
be seen by carrying out tbe di^laced. oscillator tirans — 
formation Q. = 0, 4 - (V/2511) S on tbis Haisiiltonian* ^ ^^^^ 
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H = N + I ' ( A - VQ) n. , (111.3.6) 

X ^ 

to describe tits systen. In eqn . (ill. 3.6) n^=0 relei^ to 
tbs low-spin state and n^=1 to the high-^in state, Q is a 
Erasure of the distance between the ion in question and its 
neighbotxrs, 5 is an elastic constant and (A- 70,) is the 
energy gap between the two spin states. The nodel for low- 
spin -high-spin transition described by eqn. (111,3.6) shows 
exactly the sane behaviour as the Chestnut nodsl. Thus the 
Hodel does not give any information about the change in 
symetry. It also predicts a zero high- spin population at 
TssO E. The saturation value of high- spin to low-spin popula- 
tion ratio predicted by this model for the case of cobalt (Hl) 
is 15 while from sec. III.2 we see that in cases of satura- 
tion the observed ratio is 'ul . This model gives a good 
qualitative description of the spin-state transition observed 
in transition metal c<Maplexes. Bari and Sivardiers intro- 
duced magnetic interaction between high-spin states, into 
the Hamiltonian (eqn. 111.3.6) and oaseived interesting phase 
behaviour including 'heat magnetization'. * ^ However, the 
model with magnetic interactions is also incapable of explain- 
ing some of the basic features of the transition in cobaltites 
( see section 1 ) . 

3ari,^^ Bari and Sivardiere^^ put forth a two-sublattice 
model to explain the change in symmetry observed in cobaltites. 
In the notation of eqn, (111.3*6), the model ccei side red was. 
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= N + I (A - VQ) n. + I (a + YQ)n^. (ill, 3*7) 

ieA ^ iEB 

Tais !?.odol, in addition to the phase transition between low- 
spin and high-spin states, also shows a first or second order 
phase transition to a two_sublat tics structure. This ordering 
transition is fallowed by a first or second order (if the 
preceding transition is second order, this transition is 
necessarily second order) low -spin-high -spin transition. 

This prediction is contrary to the experi rental observation, 
where initially we observe a low- spin-high-spin transition 
and this is accoinpanied or followed by an ordering transition. 


Ill .4 Models for Gpin-Stat^ Transition^ Investigated in_ the 
P_rs sen t_ Study 


111.4,1 pynapio Coupling Model. ; The models of Bari and 
Sivardiere are static models. To treat the lattice partici- 
pation in a mors cos^lste manner, ws have to incorporate the 
hinetic energy of the lattice into the Eaailt cciian of the 
model. In order to exanine the effect of the kinetic energy 
tens on the phase transition, we shall consider the following 
one- sublet t ice model: 


H = AI + B OL, 

z 

2 2 2 

^ ^ M % « 

2M+ "a" 


(111,4,1) 


B = A + VQ, 


where I is the 2x2 unit matrix, is the familiar Pauli spin 
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aatris. The eigsir states of represent the lcn-J-sp±n and high- 
spin states of the transition netal ion. Q. is the displace- 
Gsnt caused by the appropriate norrjal mode (breathing mode). 

The lo’rf-spin and higli-spin states are located at -(& + YQ.) 
and (a + VQ) respectively, , M and have their nsnal 
meaning. For the present we shall assunas that the low -spin 
and high-spin states are non-degenerate. 

The partition ftxnction of this model is given by, 

Z = Tr e“ ® (111.4.2) 

•where, from eq , (ill .4.1) we have 


^A+3 0 

E = ' K (111.4.3) 

' 0 A-3 i 


Since the operator H is diagonal , we can write the partition 
function as 


Tr 


--6(A+3) 

io 


O 1 

■3 

^-8(A-B) I 


(III . 4 . 4 ) 


= Tr . Tr,e-e<'^-®) 

Thus, to calculate the partition function, 
obtain the eigenvalues of (A+3) and (A-3) . 
■written more explicitly as 


( 111 . 4 , 5 ) 

we should first 
The opesTator (A+B) , 



2 2 
M (JJ Q 

+ - + A + Va , 


*A function of a diagonal matrix is identically equal to th© 
diagonal matrix formed by the fiaict ions of the diagonal 
elements , 
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can be trivially cast into tbe form of* hairaionic os<ri.llator 
Haciiltonian by the traaisfornjation. 



The resulting operator is 

2 ^,2 2 
M<4 Y 


(A+B) ’ = 


_ IZ 

2M 


(& - 


^ . > 
[2U^f ^ 


(III. 4. 6 ) 


( 111 . 4 . 7 ) 


Hence v;e see that the eigen values of the operator (A+B) , are 
the same as tha:t of a harmonic oscillator of fundamental 
frequency with the zero of the energy shifted by 

4- V /(2M(iJ^ ) . Similarly tie eigen values of the operator 

(A_B) are the same as that of a harmonic oscillator of funda- 
mental fi*equency with the zero of the energy shifted by, 

2 2 2 

- A - V /(2Mti^ ) . ¥c can now write the partition function Z 
of eqn. ( 111.4,2) as 0^2 

(2MaJ ^ ^ 

Z = 2 (w^^) e cosh (BA) . (III.4.8) 


This partition function being analyijic, does not shcn-j any phase 

transitions. The mean displacement Q = <Q> varies smoothly 

* 

with temperature . The difference in high-spin and low-spin 
population also varies gradually with temperaturG and is given by 


<a > = - tanh { BA ) 

z ' ' 


(111,4,9) 


This is easily calculated since the displacement can be written 
as a sum of phonon-creation and annihilation operators. 
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Tfe see frois tiais tiat the spir:-sts.ts population ratio goes to 
tmity at high tsniperaturs . The nagnetic susceptibility 
calculated from this ratio varies snoothly as $ Q-tanh 
IT we consider the high-spin state to be v-fold degenerate, the 
partition Tunction of the systen then would be given by 

Z = V Tr + Tr {III.4.10) 


This partition function is also analytic. The loean displace- 
ment of the lattice, <Q>, again shcv/s smooth variation and 
ths difference of low-spin and high-spin populations also shows a 
smooth variation. These results show that the kinetic energy 
term washes out the transition in the case of one -sublet t ice 
model. 


The above result is quits general, in that any linear 
cot^ling of the crystal field strength to the lattice, is 
incapable of leading to a phase transition. This can be 
shown by considering a model with an arbitrary linear coiapling 
of the crystal field strength to the lattice. Eamiltonian 
of such a system, in the second quantization notation, can be 


written as 


H = (h \ \ V ^ > (iii.4.n) 

k 


where and b^^ 
operators and Z 


are the phonon annihilation and creation 

* h 

is a general displacement. Any general 


For pu3Cpos6 of convenience , this displacesnent is taken to be 
the displacement of the phonon cooi*d3nate. Coxtrespanding 
lattice displacement can be arrived at by the appropriate 
transf oamaation . 
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to obtain a transfomed Eamiltozaiaa , wliicb is identical to pure 
■ononon lianiltonian vitli tl.e zero oi energy sbixtsd by a constant, 
Ix we now consider the displacensnt to correspond to tne 
dis-olacement in the two -sublattice nodel of Bari and Sivardiere 
(6^,111.3.7), we arrive at the conciasion that the kinetic 
energy tem washes away the phase transition in the case of* 
the two— sublet tice nodel also. Thus, dynanic coupling pi tne 
lattice to the crystal fisld strength, without causing a fixing 
Ox the spin states, is not the correct description of the 
systCEzs showing spin— state transitions . 


111 . 4.2 Dynanic Coup ling_J-!odel with Mixing^ oi Spin States ; 

Tfe now discuss another nodel in which the low— spin and high— spin 
stetes aix and are no longer the eigen states of* the system. 

To bring about a nixing of* the low-spin and hi^ -spin states, 
we have to couple the crystal field strength to a lattice^^ ^ 
vibration in which the transition laetal ion moves off-caatre 

*In this chapter for typographical convenience h is written 
for 
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while til© ions forning the octahedral cage reraain static. 

Sudi a sjode will caaple the high- and low -spin, states since 
the site symmetry around the transition mata.1 ion is nc longer 
octaliadral. This nixing is expected to provide a non-trivial 
iHodel in which ths higher energy high-spin state is partially 
occupied even at absolute zero. 


The cixing can he described as arising from the term. 



(HI. 4 .14) 


where cr^ is the Pauli spin matrix and Q is the displacement of 
the nomal mode in which the transition astal ion moves against 
the cage. The coupling is taken to be quadratic since the 
mixing doss not depend upon the sign of ft. The Hamiltonian 
of the model is now given by, 


H = Hq I + -^ ® CTx* 

Pq M k/ 

\ = 2M + *■ I' - » A = A and 3 = a - • ■ 


(111.4.15) 


wheare is the frequency of the lattice mode and a is the 

coupling constant in dimensionless units. At this stage we 
make the following approximatidn ; 



a m 

o 

4 


2 2 
<ft > 


o 

X 


{IIX.4. 16 ) 


This approx imat ion is necessary since a transfoa^ssatioaa that 

diagonalizes A a + 3 o should commute with S so that H 

' z , , ■ X ' ■ 'O 

is diagonalized by the same transformation. 

\ ' \ : • 
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’Ihen we seek to rind a transforaiatioii that diagonalizes H, 
we rind that a rotation arotmd Y-axis can ti^-nsfom H into 
diagonal fom. In order to obtain the actual transforaation, 
U, that diagonalizes H, we need to consider only 
Since U diagonalizes li, we can write 

2 2 

. a K 0 ) < ^ > 

U{Aa^ + Ba^) U' = Go^ , 3 = -- 

U being a rotation operator for rotation around Y-axis has the 


roin 


30 


U 


cos 4),/2 - i 0 sin , 


y 


(III. 4. 18) 


where <j) is the angle oY rotation. Substituting the r.h.s, of 
eqn. (III.4.18) in the l.h.s, of eon. (III.4.17) we have 

U (Ao^ + 30^) U“^ = (A cos <t> - 3 sin 4* ) 

+ (B cos ij. - A sin 4“ ) . (111.4,19) 

Therefore 4i is dete mined by the condition that the coeffi- 
cient of 0 ^ in son. (111,4,19) vanishes. This gives G as 


O O 1 

0 = (A^ + 


(III.4.20) 


Hence the transforned Haailtonian is giiren by 


H’ = U E V 


-1 


2 2 2 

2M 2 

2 ,2 4 

^ a 14 (ij_ pot 

+ (A + - 2r > )^ o, 


(111.4.21) 
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The parti-iicai of this systss can now be easily evaluated as 


2 =Tr e“®^ = Tr Ue”^^ U”"^ = Tre"®^' 

P a op Pi 

=Z{»^).2co^{e(A + ^ <q''> “o )^> 


(111.4.22) 


The diff srsn cs in high-spin-low-spin pcpulation,<o^> is given 


by 


<a > =Tr • a /Tr e"®^ 

21 Z 


Tn U e CT U"^ = Tr U U"^I3b tr""* 

z z 


(111.4.23) 


Froa eq , (lll,4.17)> ^^ 2 : evaluated and is given by 


U 0 TJ""^ = a cos + 0 sin i . 

Z Z X 


(111.4.24) 


Thus 

—1 - ' / 

■^cf > = Z Tr e”*^" (o cos ({ 1+0 sin ) 
z ^ z X ' 

nr-1 r X m ^ -BE' 

= Z { cos 4> Tr e '*'• * 

(IXE.4.25) 

-Be' 

3y expanding e " in. a series and using the properties of the 
Pauli aa.trices, the second tein of eon. (HI. 4 .25) can be 
shown to be proportional to the stun of trace over 0 and o . 
Since the Pauli satrices are traceless, the second tena vanishes 
The first tern can be siailarly expanded and as before, by 
using the properties of Pauli matrices, we get. 
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1 ^ 2 2 1 
<a_> = Z“ cos ♦ Tr e ° {a cosli [j{A +B 

Z 2i 


- I sinh.. [j(A^ + B^)^} 


= -Z'~^ cos i . Z (oi 2 sink [j(A^+B^)^l (111.4.26) 


Using cos $ = A/(A^+3^)^. and eqn- (III.4 .22) , -we get . 


<a > = - ' ^ ■^p- 1 tank [i{A^ + B^)^3 • 


(XEI.4.27) 


i'or a karnonic oscillator, <0, > is given by, 


<q2> 


. P- 

2m 


m 


X ;^e 


Bhm ^ 

*> .v-1 




(III. 4. 28) 


2 2 ^ 

Thus {A + 3 can be written as, 

P 9 1 w ^ P 1 

(A^+B^)"^ = A (1 + --^2° e ^^“o' 


= A# 


_ l± 


f 2 Skii) 

>- s 




(111.4.29) 


How, w© have tlie expression for <o^> a-s , 


<0 > = - - tank ( 6A$ ) 

z A$ V / 


( 111 . 4 . 30 ) 


prom eqr-.(lIT.4,30) , we see tkat = “ 'm' ’ a-t /IsO implying 

non-zero kigk-spin population at absolute zero* In tke limit 

?-»-o9 , <0 > = _ SA and tbs variation of <n > witk tei:»e 27 ature is 
’ z Z ' , ■ 

completely smootk over tke -whole t6Eperatu3?e range . 


3fkil© this uodel does not predict a phase transition, it 
does show partial occupation of the high-spin states at T=0, 
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wilicli indeed is obseived in t-^e systeais belonging to the 
categoiy of HoCoO^ . This node 1 also explains the behaviour of 
inorganic complexes, where the spin-state transiticai is not 
acconpanied by a phase transition. It appears that quadratic 
coiapling of the crystal field strength to the lattice displace- 
ment, with mixing of the two spin states is not the correct 
description of the spin-state transition in ia organic solids 
although this mecha-nisn may bs operative in causing partial 
occupation of the high-spin state at absolute zero. 

HI.4.3 Order Parameter Model ; In Section 2 of this chapter 
ws saw that the spin-state transition in rare earth cobaltites 
is accompanied by a structural transition which should be 
explained by a successful model. In order that we get a 
structural transition, we should arrive at a free energy 
expansion in terms of a suitable lattice displacement, such 
that the displacement becomes non-zero above a certain tempera- 
ture for the thermodynamically stable state. It is not possibl© 
to obtain an exact expansion of free energy in terms of the 
lattice displacement, if we consider the lattice d 5 m.amics of 
the system, as the displacement operator does not commute with 
the corresponding kinetic energy operator. However, from our • 
previous results we see that the kinetic energy term ple^rs 
important role in deteimiining the existence of a phase transition 
in the model. In order to include the effect of the kinetic 
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energy ■feena and at? trlae saisie iilae obtain an ^proxiiaate free 
energy expansion in tenas of tbe lattice displaceiaent , we bave 
eiaployed tbe patb integral netbod for compti ting the paartition 
function of our model system. In this method, the kinetic 
energy contribution to the partition function reduces to a 
constant and the partition fiznction itself is given as a simple 
integral over lattice displacement. 


To test the validity of this method, we shall first apply 
it to cases where the exact colution is Iciown. One such case 
is the one-sublattice model studied in sub-section 111.4,1, 


The partition function of this model is given by eqn , (111.4,5) 


as. 


Z = Tx* B 


,hi ih s . 


-j- © 


2 

"®^2M 




& -yq ) 


Smploying the path integral method (see eqn, II, A. 22 in the 
Appendix to this chapter) , this can be written as, 

..2 eh^. 2 


,MkT a f ^ +2ir“b) , ^ r- , 

\ * j 6 * 2 , cosii 4* /Q. } 4 dlQ 

(III..4.31) 


2irh 


when the system consists of IT such oc tab edra, in the limit of 
large N, the partition function can be approximated by. 


Z « 6 




where F(Q^) is the minimum value of the function I?{q) given by, 

o 

F(Q) = - ^ 1 In {2 cosh fl (a + VQ) 3 


. (HT.4,^) 

*The integral in eqn. (111.4.31) can be evaluated exactly. The 
results from such an integration are in agreemant with the 
exact results; that is, in both cases we do not ‘obser^ a phase 
transition and <0 2 >, and < q> calculated from both the 
approaches are in agreement. Ilawever, we rr&ke the approximation 
in the text to arrive at a free energy expansion in terms of Q . 
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The f^ee energy of* our nodel can noi:^ be -written as 

2 2 

S' = P{QJ = 1 ln{2 cosh 11(A+ 

(111.4.33) 

-where Q is given by the condition, 
n 


r- '1 

: ^ ! 

' dO ’ 

- -^=Q, 


= o 


21 


Or 


\ 

Mo). 

b 


tanb [;6(A + VQ^)II • 


(111.4.34) 


Pquation (111.4.34) haG only one solution and 'that is non-zero, 
as can bo seen graphically. This fact rules out a phase transi- 
tion in the systen. can be calculated easily and is 

given by 


tanh (BA) . (111.4.35) 

The difference of high-spin to low -spin population, <a >t 
the path integral approach is given by. 


<o > = - tanh rB(A + VQ )1 (111.4.36) 

Gos^jarison of Q. with and <a (ecn. (ill e36) 

•with (eqn. (IH.4.9)) slio'ws that in the liiait a 

o^lr approxiniation agrees -with the exact results. Also in 
agreement ■with the exact results is the fact that this approxi— 
iiKite approach does not predict a phase transition. Hence the 
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effect of the kinetic energy tera on the pkase transition is 
included in tlco approximation and at tbs same tine we have been 
successful in arriving at a free energy axpansion in terms 
of lattice displacement. 

Let us now consider the two-sub lattice case, where the 
displacements are such that octahedra on one sublattice expand 
v/hils those on the other contract . The Hamiltonian for this 
model is 

2 

% Mo) 2 

= I ^ 2M + 2 ^“i I + .(111.4.37) 

= i-)^ 

It was shown in Section 111 , 4,1 that such a model does not 
exhibit a phase transition. The exact values of <cr^> and 
are the same as those in the one sub lattice case and are given 
by eqns . ( 111 . 4 , 9 ) and ( 111 . 4 . 35 ) respectively. 


¥e can now proceed to solve this model using the path 
integral method for conputing the partition function. The 
partition function corresponding to a pair of contracting and 
expanding octahedra is given by. 





+ A + 


+ Tr e 


- 6 ( 2 ! + „ • 

+ Tr a - + Tr e 


Tff 2 2 

+ 2 ^- -i -VQ) 

■IT 2 2 

!■. a- Q 

K - ■ k - A - 1 - ’TO ) 


( 111 . 4 , 38 ) 
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goes to unity starting firoa O at Q =0. The left hanti side of 
sqn. (lII.4,4o) is a straight line "with slop© passimg 

through the origin, A second solution of eqn. (lIII,4,4o) is 
po ssible if the slope of the curve given by the right hand 
side of the equation near the origin is greater than the slope 
of the left h^d side. Thus, a non-zero Q results, in addition 
t o Q =0 , if 

VS. 

k T ^ - u- . (IIT.4,42) 

7 -T ' ' ' 

The solution also corresponds to a nininua of F{Q) (the 

second derivative of ?{Q) is greater than zero for all Q) . 

1/6, thus, arrive a.t the result that the order paraceter 
expansion of free energy in the two sub-lattice case indicates 
a phase transition in the systen, while exact and path 
integral calculations rule out phase transition in the saa© 
aodel. 

To resolve the discrepancy between the order paiTaoeter 
expansion results and the results fron exact and path integral 
calculations, we need to sxanine the order parane ter expansion 
in sone detail, ¥e assusse that the displacensnt s-t each site , 
Q., can be written as, 

= A + i <Q^> = Q and < = 0 (111.4.43) 

and write the expression fctr the partition function of the 
whole system as , 
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= / 


-6 + I j, . a.a J 

, i ^ ^ ^ j 


dQ^ 


( 111 . 4 . 44 ) 


lie re , we have introduced a snail coupling between tne displace- 
ments on sites i and j. Tliis coupling tends to cake Q to be 
tlie sane at eacb site. Using sqn . (111.4,43) in eqn. (ill. 4. 44) 
ws Iiave , 



+ AQ.AQ .) _ 

^ dQ n d A Q . ,(111.4.45) 

i ~ 


subject to the condition <AQ^> =0- Making a mean field 


approximation in eqn. (ill, 4, 45), we obtain. 


Z 



-6 yT4)(Q) 



(in. 4. 46) 


Sqn. (ill, 4. 46) in the limit of large II (using asys^totic 

. * 
approximation) yields , 


Z ss & 


(111.4.47) 


¥e see that the free energy expansion in terms of order 
parameter implicitly includes coupling between various octahedra 
in the mean field approximation, while exact and path integral 
calculations (where the integrals are calculated exactly) do 
not include such a coupling. This coupling, however, we^ it 
may be , tends to cooperatively expand and contract, the octahedra 

* 

coupling is assumed to be short-range. Neglecting' this term 
does not affect tlie essential parameters of the phase 
transition. 






term is neglsctsd since J. . is srall and the 
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thus giving rise to a phase transition. The exact solution 
discussed above corresponds to IT pairs of* octabsdra expanding 
and contracting independently and hence one doss not observe a 
phase transition. The order parameter expression Tor Tree 
energy in the one- sublattice case does not show the existence 
of* a phase transition. This ic^lies that the one -sublattice 
systsn with coupling between octahedra also does not undergo 
a phase transition . 

Tie can deternins the transition temperature by expanding 
the r.h.s. of* eqn . (ill .4. 41) around 6^ , Tor Q=Q , to obtain 
a series expansion of* Tree energy in Q . This expansion yields, 

' SL 


* l2 


(in. 4. 48) 


The t3ransition te!32peratiire is detei^ciined by equating the coef- 
2 

ficient Q to zero. This gives the transition tei 25 )erature as 

III 


___v_ 

IsM (ij . 


(ITT. 4. 49) 


The value of Q below the traxisltioii tenneraturs is given by 


Q 


*3 

“ ^ 2 


( - M„/) Ji 


3' It' 
V 




(111.4.50) 


The diTTerence in high -spin and low-spin population is given 


by <0 >, The expression Tor <0 >ia 
z z 
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siijh Q( A + VCl^)3 + siiali |j( A - ^5^)3 

< 0 _> = - ■ . V 

cosla Q(A + VQJ3+ cosh - YQ.JJ 

= - tanh (BA) . (111-4.51) 

Thus we see that although our system tindsrgoss a structural 
phase transition, the spin-state population variation is smooth 
and similar to the variation in a system without coupling 
between octahedra, ’ The susceptibility will also show a smooth 
variation with temperature. IT we compute the variatian in 
spin-state population diTTersnce on each sub lattice, we find, 

<a^> = - i {6A) (TIXA.^Z) 

sub lattice 

slaows a siaootli dependence on tbs order parameter, below T * 

c 

The two-sublattioe order paarameter model we have consi- 
dered shows a second order structural phase ti*ansition with a 
smooth variation of low- spin-high-spin population 37atio . 
Although the available experimental data in the literature 
on systems showing gradual spin-stats transitions is not 
exhaustive to enable a good comparison between theory and 
experiment, systems showing gradual spin-state transitions 
have indeed been found to exhibit anomalous Debye-T/aller factor 
changes indicating occurrence of stzTictural transition. Typical 
of these s 3 rstsms are ITdCoO^ , GdGoO^ , and the iron (ll) complex 
[i'e(papt)2^. 



138 


Systesis exhibiting spin-stats transitions az*e also 
known to taidergo thsmodyna-mically first order phase transi- 
tions at the spin-state transition 1 035)6 ra t ure . TJe have 
exasined a two-snblattice model with the spin-states coupled 
to the cube of the breathing mods displacsnent , The Hamilto- 
nian considered by us is. 


IK 

iL.' 


2 M 


MW. 






VQl’V. 


j, where Q^=(-)^Q. 


As before, employing the path integral approximation and order 
parameter expansion for free energy, we get the following 
expression for free energy: 

= fir “b - ^ E? (BA)n+ 2 % 


2 


( 111 . 4 - 5 ^) 


where , sat isfy 




r eh*' 




2. 


taiii 


p- ^ 


YQ 


'm 


+ YQ^).|(T j i.4.55) 


Sqxiation. (111*4*54) lias trwo solutions viz . , a n<m-zeTO 

solution and a Q =0 solution* TIius, tlis t'^Jo- sublatti^ systen 

n,. ■ 

with cubic c oupling also snows a -phs^s transition* To arrive 
at tli6 transition teiqpe nature and tba order of the pnasa transi 
tion, wa bave expanded tne r*]i.s. of eqn* (ttt, 4»54) in a power 

series in Q as follows; 

' ^ ' 
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V 2 


2h ~ “b ” e ^ 


„ kai 1^3 h ; 

■* / rti \~^ f\^ i ^ ^ Xl TT^ S 

Gosn (eA;3 + c 2' ^ 

T'ifT 


32 m 


^ a'* 'i£ ./ - ^ («>^ v^1 

^ .12 '"4m ' ^ ' ' 


"2M 


l& (alL) ^ 

( 111 , 4 . 56 ) 

Tile ccsndltion ior a first order pJiase transition is that the 

2 4 

coefficient of Q siionld vanish, coefficient of Q shodd be 

m * ta 

negative and the coefficient of should be positive at T=T . 

m . c 

These are satisfied when. 


,0h^.6 8^ . 6 B 2 

^4m ^ ^' 2880 ^ ^ ^ 


The transition tsnperature is given by 

3 "b^ 

^ = ,2 ■,_4 ■ 

c V h 

Like in the previous case, we can estimate 
in high-spin and low-spin population. <<j^ 
expression. 


(111.4.57) 
<a > , the difference 

55 

> is given by the 


< 0 ^> = - tanh ( B4) » (111.4.58) 

that is zero population of hi^-Spin state at T=0 a^id equal 
population of the tw o spin state s as T ■»■<«> . Tfe once again notice 
that the low ..spin and high-spin population diffei^nce over the 
whole system is independent of Q^. Tfe see froc our calculations 
that unless there is coupling of the spin states to the lattice 

■ ■ ■ ■■ ■■ ' ■ ■ . I 

displacement in even power, we will always have the above 
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es^M^ssion Tor <a ^> . I-Iowever, this difference in population 

on any one of the sublattices depends upon Q by the relation, 

. ' in, 


<0 


z 


> 

sublattice 


5 - tanh (ei) - 

(irr.4.59) 


Thus, an abrupt change in will lead to ordering of the low- 
spin and high-spin ions on different sublattices. 


In the two- sublet t ice nodel with cubic coupling, we can 
have either a first order or a second order transition. Systems 
where the spin-state transition is first order are known. 

However, thexmodynamics of such systeas have not been studied 
in detail to see whether the first order spin-state transition 
is also accompanied by a thermodynamically first order phase 
transition. The drawback of a model with pure cubic coupling 
is that it predicts the net law-spin-high-:^in population ratio 
to be independent of the order parameter. In real systems the 
net low- spin-high- spin population rationshows a jump at T . 

This feature can be incorporated into the model by adding an 
additional coupling term which couples the spin-sta.tes to the 
square of the lattice displacsmsnt . 

111.4,4 TWo Phonon Model : He saw under sub-section lU. 4, 2 
that coupling of the spin states to an optical (ion versus cage) 
mode leads to anon-zero population of the high-spin State at 
T=0,K. In this sub-section we shall make a preliminary sti:^y ©I* 
a model in which the spin states are coupled to two phonon, modes, 



viz., a 'brsatlaing aode and an optical aode . Tne brsatning siode 
to whiCii tlis spin states are coupled is the one in wiiich the 
octahedra on one sublattice expand whila those cr. the other 
contract (two-sublattice nodel) . The coupling to the optical 
mode aixes the spin states and is brought about by a quadratic 
tera in the o-i-center displacement of the transition metal 
ion. The system we shall study is described by the following 


Hasilt onian j 


= I!( 


'cC **“ • 

X X 

2K 2M 


2 ^2 „ 2.^2 

Q.- ^"^4 

^2 ■' + " 2 ) 3: + ( A + VQ^) 


X i 


(III .4 .60 ) 


= (-) q and 


X 


<ij^ and 0 )^ are the frequencies of the bz^eathing and the optical 
mode, Q is the breathing nods displacement and X, the optical 
mode displacement . 


The transformation defined by eqn. (IU.4.18) does not 
diagonalize this Eamilt onian because of the presence of the 
linear term in Q, To obtain preliminary results from the 
model of eqn . (UI .4 ,6o) , we shall the following approxi- 

5!at i ons ! ( i) we shall approsimat e t he off-diagonal term in the 
model by writing, 


Sven after making the approximation x = <x >, diagcnaiization 
by the transformation oper«.tor U of equ. (ill. 4. 18) is not - 
possible since U now becomes dependent on Q and thus does not 
commute with, H , the haimsonic o scillator part of the 
'Hamiltonian. ' • /' 
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2 

a x 


2 

a <x > a 1 

X ’ 


(III .4 .61) 


and (ii) we shall fersat "ihe node! as classics-l while consider- 
ing ths coupling to the breathing r.ode. That is to say, in 
the estination of* the partition ftinction, tjg can independently 
integrate over Q and p . Tha integration over gives a 

*4 

constant. This approx ixi^t ion can be justified since in ths 


path integral approxisation , a lin'^a.r tem in displacenent , 
in the potential function does not contr-ibuta to quantua 
mechanical correction (quantun laechanical correction is propor- 
tional to the second derivativa oi tha potential tahen with 
respect to displacement) and a quadra.tic term givss rise to 
a constant correction. Hanca wa saw in tZae linear two*-sub- 
lattice and ons-sublattice case, ths partition xunotion was 
classical except for a proportionality constant, further the 
agreement between results Trom the path integral calculation 
and axact calculation shows that classical model is a good 
approximation. It was also shown that any discrepancy between 
the order parameter expansion results and the exact results 
was due to the inherent assuc^tion of coupling between 
octahedra in ths order parameter approximation and (the 
discrepancy) does not have its origin in the path integaral 
approximation. Hence ws ars justified in approximating the 
Hamiltonian of eqn. (III.4.6 g) by 



H = I '( 2?;^ 



.. -2^2., V ■ "v: ■ 

Mta. Q. 2 I 

+ • )!+(&+ VQ^) 0^ + a <x^> 0^ ^ 


(Hi .4 .62) 





This Hamiltonian is now diagonalized by the transformation of 
eqn. (III.4.18) as in sub-section (111.4.2). 

The partition function of the above Zamiltonian in the 

case of a pair of octahedra _ one contracting and the other 

expanding - can new be written as 
' q2 

_e( P ) 1 

Z-vz {a ) 6 ^ .2 {coshg Qa + VQ)^ + a^<x^>^X 


+ cosh sf^A -VQ)^ + } dC-, 


( 111 . 4 . 63 ) 


Z(to^) is the partition ftmetioin of the optical mode. "Je see 

that the lin a-)- o, the partition runction reduces to the linear 

two-sublatt ice partition function in the path integral 

approximation ( . (III.4 .39) ) * 3n the lin the results 

are identical to the results of the modal studied in subsection 

( 111 . 4 . 2 ). Making an order parameter expansion 'wo find that 

ths free energy is given by 

2 2 

M m Qi r- 

= ^(O) + ■ 2 ^ ° ^ In {cosh SfjA + 


2 2 2^'2‘ ^ „ H- A '2 2 2 2l'5’, 

+ a <x > f + cosh g r ( A - ^ ^ < x > / , 


(HI. 4 .64) 


where Q satisfies ths equation 

El 

2 1 sirh S Wa+ 


^ “b = - 1 


L 


•V A \2 2 2 2 ^ 

VQ.^) ^ a <X > F 

m { 


2 2 21^1 


V(A - sinh |6 (a - VQ^) + a^<x 


m' 




-.1 


- f-* 2 2 2 21 ** ^ 2 2 2 21 ^ 3 

X ' coshg 1 (A + VQ ) + a <x >■ ]^ + coshes(A-VQ ) +a <x > « ^ 

i' ^ P ' ' ' M': a 12 tS,' 
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¥e see tliat Q =0 ic a solutrion of tie above eouation. For 

Q “ 

large values oif Q , the r.b.s, of the ecuatioji beiaves as 

a 

tanh (BVQ^) and hence if the slope of the r.h.s. is greater 
, 2 

than Mw, a second solution for 0 exists and this solution 
b 

■will be a non-zero solution of the equation. This indicates 
the.t the systen undergoes a phase transition. The transition 
temperature is obtained by the condition that the r.h.s. and 
l.h.s, of the expressions in sqn. (111,4,65) have the sane 
slope near the origin. This condition gives the transition 
teir^erature as 


V 




t 


(111.4,66) 


■which is the sane as that given by eqn. (IH.4.49), This 

phase transition could be second order, in -Hiiich case it agrees 

with the phase behaviour of the linear two-sublattice model 

of the previous subi? action. To see whether the system can 

show a first order phase transition, the free energy of eqn, 

(111,4.64) should be expanded in a power series in and as 

4 

before, we should see whether the coefficient of is negative 

6 ■ ■ 

while the coefficient of Q, is positive, at T . The evaluation 

■ m c 

of the coefficients is rather tedious and has not been carried 
out for the present . I-Iowsver, spin-state population difference 
can be coH^u ■feed easily, in the same manner as in subsscti<m 
(111.4.2), The expression for <a^>is given by 
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<a > = - « 
z I 


^&,+VQ^) s:kilie'_{A+VQ^)^ + a^<x^>^ 


; rj „„ ^2 2 2 2^i 

;(A+ \rQ^) + a <x > 


(A _VQ__) sinh 4 i^- VQ )^ + \ 


; I 2 ? 2 2^^ 

X [cosli 0 i(a+VQ^) + a '<x > • + 


T Jl-* 2 2 2=2 

cosli Sj (A-Va ) +a <x > » 

t_ -i-' 

(111.4.67) 


¥s Gee tliat <(3 >t3nds to zero as B-^-O or Also as T -♦■ 0 or 

z ^ 

g -♦■<», <o^>is given by 


<02.> '^- 


(A + VQ_) 

- ^ m ...^ . f .. r^- _ % -p- ^ - -w ^ - 

n x2 2 2 2l'2' 

f (a + + a <x > I 


(in . 4 . 68 ) 


This in^dies a non-zero population of the high-spin state at 
absolute zero. Another noteworthy feature in eon. (111.4,67) 
iG that the net difference between the two spin-state popula- 
tions depends on the order parameter, a fea.tare which was 
lacking in the earlier nodels. This nay give rise to more • 
interesting variations of ihe magnetic e^sceptibility with 
temperature . 

* ■ 

The results presented here afs not coa 5 >lete and the model 
itself may be stiadied in the quantum mechanical frameworfe. 
Studies in this direction are being pursued at present. 




From the Tanabe-St2g2.no diagrass, it is known that the 

4 *5 6 7 

grotind state of a transition metal ion in d , d"^, d or d* 

Gonf iguration changes from one spin state to another as the 
crystal field strength is varied. Ta solid state systecss 
where the crystal field strength is close tc such a crossover 
point, a variety of interesting m^netic and stimctuoral beha- 
viours is observed. Such observations include abznrpt or 
gradual change in the ratio of high-spin state population to 
the low-Spin stats population, a plateau in the inverse suscep- 
tibility versus temperature cuive ard change in symmetry at ths 
Spin-state transition temperature. These transitions are 
thermodynamically first or second order. Sxaoples of systems 
•undergoing ^in-state transitions are the Fe ’’’ complexes like 
lFe(phen)2(NCS)2 3 ®nd [Fe(ph6n)2{lTCSe) l^n&s , I-SriP » Co^O^ and 
the family of rare earth cobaltites, LnCoO^ (in. = rare earth 
ion) , 


Bari and Sivardiere put forward a theory to explain 
Spin-state transitions, particularly in the LaCoO^ system, ' 

Their theory does hot include the dyhanics of the lattice and 

■ ■ ■ ■ ^ 

besides does not explain the observed changes such as the ratS^ 
of spin-state populations being unity in the high temperature 
phase and the nan-zero population of the high— spin states even 
at absolute zero . In view of this , we have examined some models 

for ASpih-st ate tr^sitions in solids. ’IJe observe that 3 ard and 

; ; ^ a'*'*' h, ' , :h: ,: 
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Sivardisrs's on«- and two-sublattice node Is do not snow a phase 
transition even when a kinetic ensiigy tem is added to the 
Eani It onians, provided we do not introduce cot^) ling between the 
octahedra. The spin-state populations vary snoothly with 
teiaperatuie . The inverse susceptibility shows smooth variation 
with ten^jerature j while the behaviour at ve37y low and very 
high temperatures match the observed behaviour, the plateau 
at intermediate teiiperature is not reproduced by the model, 

¥e have examined a model in which the spin states are 
coupled to a mode in which the central transition metal ion in 
the octahedron moves ofr-center. This coupling is assumed to 
provide a mixing of the spin states since the symmetry of the 
octahedron is lost when the central ion moves off-center. This 
model does not exhibit any phase transition although there is 
some population of the high-spin state even at absolute zoix). 

To study the spin-state transition and Hie acecsapanying 
structural transition, we hscve arrived at a free energy ex^>an- 
sion in lattice displacement , This cannot be done exactly 
because of the presence of the kinetic energy term in the 
Hamiltonian. But, by using a standard path integral approxima- 
tion for evaluating the partition function, in which the kinetic 
energy term is treated approximately, we could write tl» partition 
function as an integral over displacement tally. Later by 
making an order parameter expansion of the free eneagy, in the 
limit of la»ge ®r, we have succeeded in obtaining an expan sicai 
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of* free ©nei^y ±21 "feenns of latiiice displacement. In tJiis 
approxiiaation, the two-sublattice model exhibits a secoiul order 
phase trajisiticMti, This fact is in apparent contradiction iiith 
the exact results. The discrepancy between the order parameter 
theory and the exact results is expls-ined by a closer examina- 
tion of the order parameter method. In the order paa^aaeter 
method, coupling between octahedra is tacitly assumed, in the 
mer.n field approximation, since we use a single displaceti^nt 
parameter for all octahedra. In the exact solution, on the 
other hand, displacements of each pair of octahedra are treated 
independently. Thus, we see that a linear two-sublattice model 
with coupling between octahedra exhibits a second order phase 
transition. The inverse susceptibility variation with tea^jera- 
ture is the same as in the exact model, Tl» net spin-state 
population ratio is independent of the order parameter while 
the same on each sublattice depends smoothly upon the order 
parameter. 

The two- sub lattice model with cubic cot5>ling between the 
spin states and the displacement is capable of exhibiting a 
first order phase transition. However, the net spin-state 
population ratio is again ind^endent of the order parameter. 

3 y inspection it is seen that this ratio can depeiki on the 
order parn.aeter only when the cot^ ling is to an even ptjwer in 
the displacement. 
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¥o laave investigated a two phonon aodel in which the 
coupling to off-center displacement mixes the spin states while 
the coupling to breathing mode is linear and without mixing. 

The breathing mode displacement is ccxnsidered to alternate in 
sign from site to site. This model is solved by assuming it 
to be classical with respect to the breathing mode. This is 
justified from the results on the models of subsections 
(ill. 4 .2) and (111.4.3). This model shows a spin-state transi- 
tion accompanying a structural transition. Besides it predicts 
a non-zero population of the high-spin state at absolute zero. 

The net spin- state population ratio is also dependent on the 
Older parameter - a feature which can explain the observed 
first or s.icond order spin-state transition associated with 
the structural transition and anomalous a^gnetic susceptibility 
variation with temperature. 

The models that we have considered are sia^le two and 
four state models and do not include spin -orbit coupling or 
charge transfer between different spin states. However, our 
models do show the existence of thermodynamically first or second 
order phase transitions aecos^vanying spin-state transitioiiis . Scxme 
of them also predict a non-zero population of the high-spin state 
at absolute zero and a ratio of the spin-state populations tend- 
ing to unity at high temperatures. The last model deserves to be 
examined more closely and part icularly , for variation of inverse 
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siagne tic suscsptibility -with tes^jerature. None of oiir models, 
however, ca2i predict the unusual higher temperature lower 
sysnatry phase observed in. the rare earth cobaltite systems. 
Before cornrenting any further on. the applicability and limita- 
tions of these models, it would he necessary to obtain well 
defined data on the thernxiynanic and crystal structure changes 
accon^janying spin-state transitions. Ba^jerisental sttxiies in 
this direction are now in progress in. this laboratory . 
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APPENDIX 

path INTSGRA.L method IN STATISTICAL MECHANICS 


The partition function of a systea can be ■written as the 
trace o'ver the density matrix, p , of the ^stem; 

Z = Tr / p (x,x) dx, (III.A.1) 


wlaere p(x,x) is the diagonal elenjent or the density nsatirix# The 
density Tuatrix element p(x 2 ,x^) is given by 

* -8E. 

p(x 2 ,x^) = '^i(^i) ® (III.A,2) 


■where <P^s are the eigen functions of the system and the 
corresponding eigen value. The form of the above expression is 
similar to that of a kernel for the probability of taking a 
system from (x^,t^) to viz., 

-i/h 


K(x2,t2j x^,ti) = I ® 


(III .A. 3) 


The close similarity be ■tween eqn. (ill. A. 2) and eqn. (III-A. 3 ) 
can be seen by making the substitutions u^=gh and n_j =0 in the 
earlier equation. 

The kernel of eqn. (ill .A.3) can be represented by a path 
integral of the form 

x„t,) = I .iAsg(tO 

all paths 2 
2 i/h LJ L(x,x,t)dtl 
= I e ^ Dx(t), (III.A.4> 

1 



15 ^ 


where sjx(t)3 action for a path connecting (^ 2 *^ 2 ^ 

{3c^,t^), L{x,:i,t) is the Lagrangian of the systen and D3c(t) inplisi 

integration over all paths connecting * Since 

^(x^jt^; x^,t^) and p(x 2 ,x^) are similar in form and also can he 

* 

snown to satisfy similar equations of motion , 
represented as a path integral by performing the change of 
variable from t to -iu, in ecn. (ill, A. 4). Thus the density 


matrix element can be written as 

h 


= fe 




Dx(u) , 


(III .A. 5) 


where x^ is the differential of x w.r.t U. The path integration 
is performed over all paths connecting the points (x^,0) arid 

, 6h). For systems where Lagrangian can be written in ihe form 

L = I x^ - V(x) , (III.A.6) 

the change of variable from t to -iu gives the 'Lagrangian' of 
eqn. (ill. A, 5) as 


u 


m 

2 


ij" - n^) 


(III-A.7) 


Thiffi the density matrix element of such a system can be written 
as 


p(x2,x^)=Je^ ° ^ ^ 


(IU A. 8) 


The partition function now becomes 


In the case of the density matrix, -i&h assumes the role of 
time , 
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^1 ^ C J (§ x^+v(x))dun 

Z = ;p(x^,x^) dx^ = Jdx^ / g ^2 u ^ Dx(u)l 

^ (III. A. 9) 

Tlie paths over which integration is perforrisd begin and end at x^ , 
At high teaipeiratures , the 'tine’ {6h) in which a path has to be 
coiapleted is snail. Hence longer paths will have higher 'kinetic 
energy’ . Consequently they do not contribute significantly to 
the integral. Since paths which contribute significantly are 
shorter, the potential V{x) can be considered t o be a consteint 
for the entire path. The path integral that need be evaluated 
now is just a free particle path integral and can be evaluated 
easily. The resulting partition function has tbs forai 


= ( 2) J® • 


(III .A. 10} 


For our purpose, we need to estiaate the partition function 
over all temperatures. At low temperatures longer paths do make 
significant contributions to the path integral. So, we no longer 
can treat the potential as a constant and should take into 
account the change in the potential along a path. ¥© can approxima- 
tely treat the variation of the potential along a path, by 
expanding the potential around a point x in Taylor series. For 
convenience and accuracy, we can take this point to be tl» mean 
position of a path, x, defined by 

1 6h 

X = / x(u) du. (lII.A.It) 

o 

The path integration is carried oat over all paths which have the 
sane x. A second integration over x will give the pari; it ion 
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function. Expanding the potential in Taylor series upto the 3 rd 

tenn, we can write 

gh gh 

J V g(u)3iu = gh V(x) + I g(u)-x3v’(x) du 


+ j j [x(«) - 5 II (III.A.12) 

o 

Ths second term in the expansion vanishes from the definition of 
X (eqn* III. A. 1 1 ) . Hence, we have for the partition function 





i j (f 


.2 

r + 


(x(u)_x)^ -■ _ 

O v‘'(£)}duf 1 

^Dx(u)| 


(III .A. 13 ) 


Y(x) is a constant since x is a constant for the set of paths 
considered in the path integration. 

The constraint, eqn. (ill. A. 11), can be incoirporated into 
the path integral by rewriting it as 


gh 


4 (^-^) = 


0 


(III.A.14) 


Substituting y = x - x gives the constraint in the form 
gh 

3 y du = 0 . 

o 

The path integral then as sun® s the form 
x. -x f B 1 


(III .A. 15 ) 


P.I. = j e ^ -^Dy(u) . 

X —X 

“ (III.A. 16 ) 

Tho constraint, now written in the form of eqn. (III.A,15) cto 
be incorporated into the path integral in the following way. 


y^ ‘ill! 
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Since all the paths in the path integral have the saae average 

eh 

position, y=sO, we snultiply the path integral by 6( / y du) . Now, 

O ’ 

we use the integral representation of the 5 jPunction viz,, 

.eh 


eh 


6(/ y du) = J 


„ ik ^ y dn 
® 2ir 


(III.A.17) 


Thus the path integral with the constrair.t noi-j becosnes 


P.I. v;ith constraint 


= J 


dk f 
2t ■' 




e 


ir- 31 .2 

"h' o ^2 ^ ■*‘2 


-X 


^ y^ ^^"{o) - ihky}du i 

Dy(u) (III.A.18) 


The path integration in eqr. . (ill. A. 18) is periorried over all 
paths beginning and ending at -x. To complete the path 
integration, integration over x^ (to include all initial points) 


also has to be carried out. 


The path integration in eqn. (ill .A. 18) is carried out by 
first writing the paths y(u) as 

y(u) = y^ 3 ^(u) + g(u) , (iil.A,i9) 

where y is the classical path corresponding to the Lagrai^ian 
Ox 

±n thjB exponent and g(u) gives the diiiex'ence between the classical 

path end the path in question. The end points of g(u) are (o,0) 

as y satisfies the required bo-undary c audit ions. Substituting 
OX 

eqn. (ill .A. 19) in eqn. (ill .A. 18) , the latter equation can be re- 


written as 
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p.i; = [ I 


^ dx, e 
2tr “^1 


-I I ° [ni 


(III.A. 20 ) 


where S is the classical action. First order tazrss in g and g 
Ox 

do not appear in tlie exponent oi* tiie abo^ve equation a,s tlie 

classical path is deterTdned by setting the first oixiez* correction 

to zero. The path integral over g can be easily evaluated since 

g can be extended odd periodically in u (this is possible since 

g(u)=:0 at u=0 and u= 6h) *with a period 26li, The function, g can 

then be expressed as a Fourier sine series. Integration over all 

the paths, g{u) , is novj reduced to a set cf simple integrals 

over the Fourier coefficients ani the intagrations yield a 

1 

constant times (6h)^, Integrations over k and x^ can be performed 

after expanding S correct to first order in The path 

bJL 


integral is finally given by 


P.I. =/' 


2 irh sin (tuHi) 




o2, 2 

- ^ . V» (x^ 

• rtOcT A 24m 

. p) ® 


(III.A.21) 


for snjall V’ (x) . Substituting this in eqn. (ill. A. 13 ) gives the 


partition function as 


(m_)t I 




(HI, A. 22) 


Thus, we see that in the expression for the partition function, 
only the distance X appears and the kinetic eneigy contribution 
to the partition function appears as a correction to the potential. 
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